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THE  PROBLEM 


Deter  nunc  .  md  study  lie  aeouslic  pi  ope  i  lies  ot  sea-tlooi  sediments  and  acoustic 
models  ot  the  si* a  tluui .specifically.  use  measured  values  ot  comprcssional-wavc  (sound) 
vol* u'll y  and  density .  will:  computed  values  ol  sednnenl  hulk  modulus  to  compute  olhei 
elastic  constants,  and  do  ter  nuiio  appropriate  elastic  models  lot  the  sea  lluoi 


RESULTS 


! .  Marine  sediments  can  be  considered  as  elastic  media  and  the  equations  ot 
llouheun  elasticity  cjn  be  used  to  compute  unmeasured  elastic  constants.  When  sound 
attenuation  must  he  consideied,  a  viscoelastic  nu:del  is  I  a  voted  in  winch  the  !  i  *.  i  O  1 1  \ 
modulus.  p.  ami  Ijiiic's  constjnl .  A.  m  the  equations  ol  elasiiciu  ;ko  replaced  b>  coin- 
ple.\  I  .ante  const.mtstp  a  ip'  ami  A  +  i  A  ’ )  winch  are  independent  ot  t  requeue)  In  tins 
model,  p  and  A  represent  elastic  response,  and  ip’  and  rA'  represent  dumping  ot  wave 
eneigv  I  his  model  implies  that  11.110  velocities  and  the  ilissipation  luncium.  I  (/.  are 
independent  ol  frequency.  and  attenuation  in  uli  pet  unit  length  varies  linearly  with 
fiequeiKs  in  the  range  from  a  lew  11/  to  the  megahertz  range 

2.  Density  jttd  eompressioital  wave  velocit)  were  measured  in  the  present  studs' 

A  computed  value  lor  the  system  hulk  modulus  (follow mg  (.iassmann.  1051 )  w  as  nseJ 
as  the  thud  constant  required  to  compute  the  other  elastic  constants  with,  the  equations 
ol  elasticity  .  'I  hese  computations  wcie  thus  based  on  llieoiy  without  empmcal  factors 
01  constants. 

The  components  ot  the  computed  system  hulk  modulus  are  porosity  .  the  bulk 
modulus  ol  poie  water,  an  aggregate  hulk  modulus  ol  mineral  grams,  and  a  hulk  modulus 
ol  the  structure,  m  Iramc.  tormed  by  1  lie  mineral  grams  (mod  values  lor  the  hulk 
moduli  of  distilled  and  seawater,  and  most  of  the  common  minerals  01  sediments,  have 
been  established  in  recent  seats,  this  leaves  only  a  value  for  the  frame-bulk  modulus  to 
compute  a  hulk  modulus  for  the  a ater- mineral  system  ol  the  sediment 

Curves  and  regression  equations  relating  porosity  and  dy  namic  frame  bulk 
modulus  were  derived  and  used  in  computations  ot  a  system  bulk  modulus,  which,  with 
measured  and  vclcc*!'  \«  js  U'cJ  *  ■„*  compute  other 

0 1 i* s r ic  constants 

v  The  computations  ol  elastic  constants  as  discussed  in  tins  report,  compared  witli 
other  laboratory  anj  i>i  stfu  measurements  (such  as  shear-wave  velocity )  indicate  that  liic 
equations  o!  elasticity  can  be  used  *o  derive  reasonable  values  tor  those  elastic  constants 
not  measured,  It  so,  the  listed  values  m  (Ik*  tables,  and  the  method,  can  be  used  to  pre¬ 
dict  these  constants  tor  the  major  rcdimcm  types 
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sediments  l  lie  aU'Mm.ih"ii  «*l  vonipiesMen.il  (  Mum*i  J  ’ajux  1  he  velocities  o!  shc.u  u.incs 
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\.u i.it ions  \--t  propeilies  vmiIi  depth  in  the  sea  limn. 
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Mippoit  undewaki  sound  pi up.ip.it  ion  and  i*coph\  sie.i!  studies 
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PREFACE 


IliiMepott  iNl’jit  I!  m  j  m'ik'miI  lccliiiiv.il  Publications  on  sound  velocity ■. 
cIjmiciI)  ,  ami  i  chi  led  piopeitiesol  tn.nino  sediments  lioin  (luce  major  etiviionn.enis 
»'l  1 1 1 C  N'Ulll  r.KllK"  Ilk*  OMlIl1KMI1.l1  lCI  I  .ICC  (  shell  .nld  nI  ope  ).  a!  »\  \Sj|  phi  111  ( til  i  bid  )  I  0  ). 
am!  ahs  nn.i!  lull  (jk'lapic) 

Pail  1(1!’  I -1  a )  detailed  the  incaMiicnurni  and  computation  o!  mass  p! i n  sical 
pi  open  ion  oi  the  sediments,  inclinum*  -tensity  .  \  H%»v*  1 1  >  .  pi  ain  size.  impedance,  icllvction 
coellic'.oiiis.  and  hoi  tom  loss 

1  III’  piCsClll  K'piM  l  (IP  II  1  J  pi  Cm' I  It  N  ili'VUSM'  »lis  o!  cLinI  'C  .tiki  \  INCOChlsI  1C 
models  Ini  v  aler  salinated  pm  oils  media,  and  meastii cinents  ami  computations  ol 
elastic  constants. 

Pait  ill  (  I  P  1-Js>  will  he  coiiicrnci!  with  predictions  and  vonipiit.it ions  .»l 
in  \tfu  physical  propet  tie v 

1  liioujdiotit  cavh  icpmt,  lelercnces  aie  nude  to  *lie»»thei  two  studies, 
as  appiopiMie. 


REVERSL  SIDE  BLANK 
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ROCKS  AND  SEDIMENTS  AS  ELASTIC  MEDIA 
Introduction 

M  .1  ■  i  *  p.ijVl'Mi  lliO  il'CM-IIC  “I  ^Ot'I'^V  .  00>»p!:\  -ICs.  pb\  Ml',  .rill  mmI  MIC .  I  i  .i  I  i  !i  s 
uic  il'icclci!  I-.*  1 1  iv  ijuo'lj-  i»l  vs  I  ic  I  i :  ci  walci-saitii  .nod  look  v  vim  is.  ,n*.l  cij>  '  ;iio 
ci.is'K  1'Hil‘Os  4 1 1  •  I  ci'.'-’-IC.ll  Sl'MsOiM  \  !sc*0l,lsl  li  .  <t|  «‘l  Ik'l  .  IliOsll.i  .  .itlii  'A  l-.-.'l  liCl  .  li  ll'CV  ,MC 
\  I  s-J  "0 1 .1  s  i  K  «  I  Ik'  C  IllClIl.l  .11  C  i  If  SC(  lit  I  I  -  ifi.C.ll!.  Kolv  111  \  ;  ‘Ip  i  .  M.lWSC'l.  "i  * ‘ill  01  111-'  k  I ' 

I  >151  lilt!  lliC  p.l'I  I  MJ' s  .1  VIC.l'  tlc.il  c|  1  jiC  '!  0 1  !C  1 !  .il  :-.l  0\ pi'f  SlSiCI  .1  .li  V.fD.  (Ml  •“-!!>  '  Oi 
!■  'I  Ks  i  ll.'.s  lie. ill  Will  1 1 1 1  s  1 1 1  JSY  'Cl.llOti  I}li0sl  S  ’liv  III  ill!'  J'CI  |.‘ii.  Ospci  l.tliv  SIliCl  .1 1  ■  ’  1 1 1 

II  Mu !.  c.iu'f  ;i!  ov  per  ii'it’iii  .il  \»  .'i  I  will*  '  ii  c!-’-.  •  >  •  ■;»  .  1  :ic.: i :p:;;c;:5 .  5.  -^c !  !•.  , 

will  i  1  home!  n  .il  •  I  nil  10  s.  Ii.»s  Osl;ibhv!:otl  |  lia?  i » *  i  ks  (  nit  luiiiih:  w  .ilci  s.i!  U'.iU'tl  n\  k'  J  .no. 
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m.icff  scfpii  .ilk  isitiiopio.  1  oi /.iplii  ( 1  jilopioil  il.is  p'niii  Hi  siow  in  lnsilovohip 
moil!  i*l  l lie  lhom\  ol  chiishIivI.iIihu.  Iii  !‘MI.  Ui«»l  wp'lc.i  cl.issu  p.ipci  fn  thioo- 
■JiinonsioiKil  ooiisfhtl.il mil  m  winch  lie  adopted  i lie  s.mho  \iowp hi:',  llns  p.ipci  lias  imi 
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received  the  attention  it  deserves  Hi.it  a. id  Willis  f  |957)*  showed  how  the  elastic  con¬ 
st. nits,  denved  in  I 'Ml .  could  he  experimentally  determined  by  static  tests  on  jacketed 
and  iiniucheied  samples  ot  poious  media.  In  |9So,  hint  extended  ills  analyses  to  coin- 
pressional  and  slie.n  waves  by  adding  menial  terms  to  Ins  theory  ol  three-dimensional 
consolidation,  he  later  extended  his  concepts  to  viscoelastic  and  anisotropic  bodies.  In  a 
recent  papet  (Riot,  |9(,2)  he  sinnin.n ized  and  extended  much  ol  Ins  work  ot  the  previous 
JO  y  ears 

In  Ins  various  papers.  Hint  discussed  the  lull  range  of  sy  stems  in  which  water 
within  pore  spaces  does  or  does  not  move  with  the  solids  upon  imposition  ol  a  small 
stress  such  as  that  of  a  sound  wave.  In  some  acoustic  models  and  theoretical  siudies,  this 
movement  or  flow  of  water  through  the  sediment  mineral  structure  has  beer,  considered 
to  be  of  the  I’oisuclle  type.  In  the  last  several  decades  it  has  been  determined  that  the 
simple  flow  equations  of  the  IVisuelle  ty  pe  (derived  from  flow  of  water  through  tubes) 
do  not  hold  for  real,  in  situ  sediments.  These  equations  have  to  be  considerably  altered, 
even  lor  clean  sands,  and  are  not  applicable  to  relatively  impermeable  clays  (Yong  and 
Warkentin,  1 906)  The  question  of  relative  water  movement  is  a  critical  key  to  whether 
or  not  the  equations  of  elasticity  (f  e. ,  llookeatt  respi  nses)  can  he  used  in  studies  of 
rocks  and  sediments.  If  the  pore  water  does  not  move  significantly  with  respect  to  the 
solids,  then  the  effective  density  of  the  medium  is  the  sum  of  the  mass  of  the  w  ater  and 
solids  in  a  unit  volume,  water  viscosity  need  not  be  considered,  and  the  equations  of  ideal 
elasticity  can  be  used  This  is  the  "closed  system"  of  Gassmann  ( 1 95  I ) 

To  avoid  an  extensive  review  of  the  elastic  relationships  in  porous  meuia  (sedi¬ 
ments  and  rocks),  the  following  concepts  (believed  by  the  author  to  be  true)  are  slated, 
with  references  which  support  or  have  an  important  (affirmative)  bearing  on  the  concept 

1  If  a  saturated  porous  medium  is  a  “closed  system"  it  resf  ids  to  small  stresses 
(such  as  that  of  a  sound  wave)  as  a  llookeap  body  (strain  :R  recoverable  upon  removal  of 
stress)  Irt  a  linearly  elastic  sy  stem:  any  deviations  front  llookeatt  response  are  so  small 
that  they  can  he  neglected  for  the  purposes  of  undei  water  acoustics  and  geophysics,  and 
the  sy  stem  can  be  studied  with  the  equations  of  classical  elasticity  I  btol.  1941,  1956. 
1962;  Biot  and  Willis,  1957.  (iassmann,  1951  ;  Morse.  1952;  White  and  Sengbush,  1953; 
Zwikker  and  kosten,  1949. Brandt.  1955;  Hamilton  vial.,  1956;  Paterson.  1956;  Wy  Hie . 
Gregory  ,  and  Gardner,  1956;  Laughton,  1957.  Jones,  1958;  Knopoff  and  MacDonald, 
I95H;  Me1  tonal  ct  al  .  1958.  Birch.  I960.  1961,  1966:  Barkan,  1962; Geertsma  and  Sinit. 
|96l  ;Nafeand  Drake,  I9(>3.  Brutsaert.  1964;  Simmons  and  Brace,  1965.  White,  1965). 

2.  Must  rocks,  composed  of  randomly -oriented  mineral  crystals,  are  macro  cop- 
ically  isotropic  to  congressional  and  shear  waves  in  a  frequency  range  from  a  few  Hz  to 
the  megahertz  range  if  the  grain  sizes  arc  much  smaller  than  wavelengths;  the  elastic 
constants  "t  such  systems  can  be  computed  with  the  Vojgt-Reuss-llil!  avei aging  method 
(Hill,  1963),  discussed  below.  A  recent  excellent  summary  oT  the  evidence  for  the  Voigt- 
Reuss-llill  averaging  method  lias  been  published  by  Anderson  and  Ltebermann  (1968) 

In  addition  to  references  cited  in  tiie  previous  paragraph,  attention  is  called  to  reports  by 
Brace  ( 1965 ).  Peselnick  ( 1962),  and  Christensen  ( 1965.  1 966) 

3.  Static  measurements  ol  the  elastic  constants  ol  a  saturated,  porous  rock  or  sedi 
ment  arc  essentially  isothermal:  dynamic  constants  are  adiabatic  In  rocivS  and  minerals, 
isothermal  and  adiabatic  measurements  can  be  experimentally  and  theoretically  related 
in  a  linearly  elastic  system  (Biot  and  Willis,  1957,  Biot,  1962:  Anderson  ct  al.,  1968, 
Brutsaert,  1964,  Trietei,  1959,  Anderson  and  Schrciber,  1965,  and  othcis).  In  connection 


' Sc v  lid  of  references  at  end  ot  report 


8 


with  tins  concept  ami  tor  discussions  of  the  tael  that  elastic  constants  measured  by  static 
methods  are  smaller  than  those  measured  by  dynamic  methods,  attention  is  directed  to 
reports  by  Simmons  and  Brace  (  l%5 )  and  Walsh  ( 1965) 

4  Liquid  salination  ol  porou  rocks  and  sediments  increases  the  velocity  ol  com 
press lonal  waves,  and  decreases  the  velocity  ol  shear  waves  ((>assmann,  1951;  W  lute  and 
Sengbush,  l *-#5 3 ;  Utol ,  1956, King,  19ot>)  Ligwrcs  12  and  I  .»  illustrate  these  phenomena 
(Shell  Development  Co.,  private  coiiiiiinnicjlioii,  1 90S). 

5  It  is  not  neeessai  to  postulate  movement  el  viscous  pore-water  relative  to 
mineral  grams  to  explain  attenuation  m  satmaieJ  porous  media.  Olitei  sourd-.'iiennaiinn 
1  actors  not  involving  sediment-structural  permeability,  water  viscosity  and  lime  aie 
apparently  present  for  example,  thermal  losses  and  losses  due  to  internal  Inetion  tsco 
Bradley  and  I  on.  19(>6.  for  a  leeeni  resume) 

t>.  The  velocities  of  compressional  and  slreat  waves  are  imlepc  \1  nt  of  Requeues 
(fe.,  no  dispersion  or  change  of  velocity  with  licquency)  funn  a  lew  I  *rtz  into  the  mega- 
hertz  range,  or  dispersion  is  negligibly  small,  experimental  evidence  on  this  subject  is 
discussed  in  a  later  section. 


Elastic  and  Viscoelastic  Models  for  Rocks 
and  Sediments 

In  the  above  discussion  of  theoretical  development,  it  was  stateJ  oi  assumed  that 
the  numbered  statements  were  true  only  for  small  stresses  In  the  field  of  soil  mechanics, 
especially ,  large  static  or  dynamic  stresses  have  to  be  considered,  and  over  the  full  range  ol 
stresses,  sediments  are  both  elastic  and  viscoelastic,  't  one  and  Warkentin  (  l9oo.  p  50-94) 
have  a  good  discussion  of  the  various  models  and  elements  within  llie  models  which 
describe  this  behavior.  These  writers  favored  the  Burger  model  for  mechanical  simulation 
ol  actual  soil  behavior. 

The  Burger  model  is  lormed  by  coupling  the  Kelvm-Voigt  model  in  series  with 
the  Maxwell  model  ('l  ong  and  Warkentin.  196b.  p  9U-93,  figs.  4,  12-15).  In  the  Burger 
model  very  small  stresses  arc  apt  to  he  recoverable  (ie  .  Hookean),  but  if  the  magnitude 
of  a  stress  is  sufficient,  the  medium  behaves  as  a  viscoelastic  body.  In  the  fields  ol  soil 
mechanics  a  :d  foundation  engineering,  the  Hookean  mode!  and  equations  are  commonly 
used  lot  derivations  of  dy  namic  elastic  constants  and  studies  of  vibrating  lojds  (e  g. . 
Barkan,  |9b2:  Heukelonv  1961.  Jones.  195S.  l-vison.  1956.  Hardin  and  Richati.  196a; 
liaii  and  Kichan,  196.5;  Kichart  and  Whitman,  1967). 

In  the  fie'ds  of  physics  and  geophysics,  studies  of  the  elasticity  of  minerals  and 
rocks  have  uemonstraled  that  the  elastic  equations  of  the  Hookean  system  adequately 
define  the  velocities  of  compressional  and  shear  waves,  these  equations  arc  conveniently 
interrelated  in  a  table  by  Birch  (  1961  .  p.  2206).  This  field  has  been  summarized  by 
Birch  (  1966)  and  by  Anderson  and  Lie1  erinann  (|9oh.(;  papers  of  special  interest  are  by 
Christensen  I  1966a.  b).  Brace  (  1965),  and  Simmons  and  Brace  (  1965) 

Although  the  elastic  equations  of  the  Hookean  model  adequately  account  for 
wave  velocities  in  earth  materials,  they  do  not  provide  for  wave-energy  losses  (attenuation ) 
m  these  media  To  account  lor  attenuation,  vaiious  viscoelastic  or  "neat -elastic"  models 
and  equations  have  been  proposed.  One  of  these,  which  has  been  studied  in  connection 
with  rocks  and  sediments,  is  the  Kelvm-Voigt  model,  in  which,  as  originally  defined, 
compressional  wave  velocity  varies  with  frequency;  and  attenuation,  at  frequencies  of 
interest  in  underwater  acoustics  and  geophysics,  increases  with  the  square  of  frequency 
White  ( i965,  p.  110-112)  has  a  thorough  discussion  of  theory  jnd  experimental  evidence 
on  this  subject,  and  concludes  I  p.  I  12)  that  neither  velocu  nor  attenuation  shows  this 


trequency  dependence,  and  the  Voigt  solid  cannot  he  considered  an  adequate  model  of 
earth  materials.  The  cxpcuinenta!  evidence  of  this  leport  (and  ot  reports  in  preparation) 
is  in  accoul  with  tins  conclusion. 

A  growing  bods  ol  experimental  evidence  in  the  field  and  laboratory  has  tenta¬ 
tively  placed  several  piobable  test! ictions  on  nonelastic  models  and  equations.  I  or 
example,  velocities  and  energy  damping  are  independent  of  frequency,  and  spatial  attenu¬ 
ation  increases  approximately  linearly  with  tiequency  ill  the  ranee  ol  trequencies  ol  most 
mteiosl  in  engineering,  inaime  geophysics,  and  underwater  acoustics  (a  lew  liorl/.  to 
the  megahertz  range), 

A  viscoelastic  or  "near-elastic  model  within  the  limits  of  the  above  restrictions, 
without  specification  ol  the  mechanics  ot  wave  attenuation,  has  been  discussed  by  White 
(  1 96s,  p.  79-154).  Bucket  (  1964),  Hamilton  ft  al  (  1969).  and  otheis.  In  this  model  the 
Lame  constants,  p  and  X.  in  the  equations  of  elasticity  are  replaced  by  complex  Lame 
constants  (p  +  tp  I  and  (X  +  i\')  in  which  p  and  X  represent  elastic  response,  and  p  and 
X'  represent  damping  of  wave  energy.  This  approach  is  based  on  the  assumption  that  all 
elements  of  the  complex  i  ante  constants  are  independent  of  frequency  in  the  range  of 
frequencies  of  interest  in  most  problems  of  engineering,  geophysics,  and  underwater 
acoustics.  The  basic  equations  of  this  and  associated  models  are  presented  in  papers  by 
[■wing  cr  al.  t  1 957 ),  Bucker  ( 1  964),  White  ( 1 965 ),  and  I  lamilton  ft  al.  ( 1969) 

The  assumption  that  the  complex  Lame  constants  are  independent  ol  frequency 
from  a  lew  hertz  to  the  megahertz  range  implies  that  their  elastic,  or  real,  parts.  P  and  X, 
are  independent  of  frequency.  Because  these  constants,  and  density,  govern  wave 
velocities  it  is  further  implied  that  both  compressional-  and  shear-wave  velocities  ate 
independent  of  frequency  (in  literature:  "no  dispersion  of  velocity  with  frequency "  ). 

The  subject  of  possible  dispersion  will  be  discussed  in  a  later  section. 

If  the  complex  Larue  constants  ate  independent  of  frequency ,  their  imaginary 
parts,  which  govern  energy  losses  of  wave  propagating  through  a  medium,  must  he 
independent  of  frequency .  Biot  ( 1962,  p  14981  noted  that  the  imaginary  part  of  com¬ 
plex  rigidity  (and  other  operators  involved  in  energy  losses)  represent  a  damping  which 
varies  very  little  within  a  relatively  large  range  of  frequency.  The  independence  of  the 
dissipative  parts  of  the  complex  Lame  constants  from  frequency  is  implied  in  the  studies 
of  energy  losses  in  earth  materials,  where  the  specific  dissipation  function,  I  ;Q.  and  the 
logarithmic  decrement,  A,  have  been  shown  to  be  independent  ol'  frequency  over  a  range 
of  at  least  I0S  Hz  (knopoff  and  Macdonald,  1958;  White,  ll,65  :  Bradley  and  L'ort,  1966, 
Allwell  and  Rantana,  1966) 

Impoilani  evidence  that  the  complex  Lame  constants  uie  independent  of  ire 
queue)  in  high-porosity  sediments  was  tecently  furnished  by  Cohen  ( 1968).  lie  meas¬ 
ured  both  g(  and  in'  in  complex  rigidity  (p  +  in  )  in  artificial,  laboratory  sediments 
composed  ol  kaolinite  and  bentonite  in  distilled  water,  with  and  without  a  defloceulating 
agent.  Cohen  demonstrated  that  both  p  and  in'  were  independent  of  frequency  in  the 
range  8.6  to  43.2  kHz  in  flocculated  clay,  but  when  a  defloceulating  agent  was  added, 
the  flocculated  structure  of  the  clay  sediment  dispersed,  the  material  lost  all  rigidity,  p, 
and  behaved  as  a  Newtonian  viscous  fluid  in  which  in'  was  “viscosity and  was  linearly 
dependent  on  frequency.  The  addition  of  35.5  ppt  of  NaCl  caused  reflocculation,  and 
complex  rigidity  w  as  the  same  as  before.  Because  in'  governs  energy  dissipation,  this 
demonstration  that  in  was  independent  of  frequency  extends  the  range  of  frequency- 
independence  of  the  dissipation  constant,  1(9,  from  rock  to  flocculated  clay  structures. 
Other  aspects  ul  this  study  will  be  discussed  m  later  sections 

Another  implication  of  the  assumption  that  the  complex  Lame  constants  are 
independent  of  frequency  is  that  attenuation  in  dB/unit  length,  a,  increases  linearly 
with  frequency, /(White.  1965,  p.  98).  Recent  summaries  of  work  in  this  field  indicate 
that,  for  most  earth  materials,  there  is  a  small  variation  around  linearity  in  the  range  of 


frequent. k‘5  ol  interest  in  underwater  acoustics  and  geophysics;  that  is,  in  the  relationship, 
a  =  k  /';  the  exponent  ";i"  is  about  one  In  clean  sands,  however,  there  is  some  evidence 
that  the  exponent  may  he  nearer  0  5  (Nolle,  er  a/.,  1963;  Hampton,  1967),  but  this  matter 
requires  further  study  in  natural  marine  sands.  A  linear  relation  of  attenuation  with 
frequency  lias  been  successfully  used  by  C  ole  ( 1965 )  in  cxtrapola'ing  resonant-chamber 
attenuations  :n  mat  me  sediments  at  20  to  40  kHz  down  to  100  Hz;  and  Wood  and 
Weston  ( 1964).  m  an  experiment  m  a  tidal  mud  llat,  found  an  approximate  linear  depend¬ 
ence  in  the  frequency  range  from  4  to  50  kHz. 

The  elements  of  the  complex  l  ame  constants  (p, /a',  A,  and  A')  can  be  computed 
if  the  velocity  and  attenuation  of  the  comprcssional  wave,  velocity  of  (tie  shear  wave,  and 
density  arc  known,  and  bulk  viscosity  is  assumed  to  be  zero  "Bulk  viscosity "  is  a 
phenomenon  related  to  the  time  lag  which  may  occur  between  pressure  and  volume 
change  (Jaeger,  1962,  p.  71 ;  Biot,  1962,  p  1492);  it  is  analogous  to  the  elastic  bulk 
modulus,  k  -  A  +  2/j : i.  belting  bulk  viscosity  equal  to  zero  (A'  -  -2/u',3)  is  a  simplifica¬ 
tion  which  has  been  used  by  Jaeger  ( 1962,  p.  71 ).  T  vving,  ci  al.  ( 1957.  p.  273).  Biot 
( 1 962).  and  others. 

Recently,  in  situ  measurements  weie  made  (from  a  research  submersible)  in 
sands  and  silt-clays  ol  the  velocities  and  attenuations  ol  compressions!  waves,  and  veloc¬ 
ities  of  Stonelcy  waves  (from  which  shear  waves  were  computed).  These  measurements 
alloweJ  tentative  evaluations  ol  the  above  viscoelastic  model,  and  computations  ol  both 
elastic  (Hookeun)  constants  and  the  complex  Lame  constants  of  the  viscoelastic  model 
(Hamilton  ct  al..  1969).  These  computations  affirmed  that  either  Hookean-clastic  or 
viscoelastic  equations  can  be  used  to  derive  the  same  eomprcssional-  and  shear-wave 
velocities  and  associated  clastic  constants. 

The  author  believes  that  the  evidence  of  the  present  report,  against  the  theo¬ 
retical  and  experimental  bachgiound  outlined  above,  strongly  supports  the  use  of  the 
Hookeun  ( “closed  sy  stem”)  elastic  model  and  its  equations  in  computations  and  predic¬ 
tions  ol  the  velocities  of  eomprcssional  and  shear  waves,  and  other  elastic  constants, 
in  iiatuial  marine  sediments. 


Pertinent  Equations  of  Elasticity 

II  the  llookean  model,  or  the  Hookean  component  in  more  complex  models, 
is  approximately  correct  for  studs  ing  the  velocities  of  compressmnal  and  shear  waves  m 
saturated  sediments,  and  classical  equations  oi  elasticity  can  be  used  to  compute  those 
constants  not  measured,  density  ,  and  any  two  of  the  other  elastic  constants  ate  required 
In  the  case  of  water-saturated  sediments,  the  density  and  compressional-wave  velocity 
can  be  easily  measured  or  can  be  closely  predicted  (as  discussed  in  Parts  I  and  111 ).  The 
best  third  constant  for  computations  of  the  elastic  constants  of  marine  sediments  is  the 
velocity  of  shear  waves,  at  least  for  purposes  of  underwater  sound  and  geophysics. 

Information  on  the  velocity  of  shear  waves  in  marine  sediments  is  rare,  in  fact, 
it  is  so  rare  that  one  of  the  principal  contributions  of  this  report  could  be  in  the  computa¬ 
tions  of  shear-wave  velocities  I  which,  in  effect,  constitute  predictions  of  the  values 
expectable  in  natural  marine  sediments).  Comparison  of  experimental  with  computed 
values  of  shear-wave  velocities  w  ill  be  made  in  a  later  section. 

Lacking  sufficient  information  on  shear  wave  velocities,  the  third  constant  used 
here  for  computations  of  elastic  constants  is  the  bulk  modulus,  or  incompressibility  ( the 
reciprocal  of  compressibility  ).  This  constant  was  selected  because  it  is  possible  to  com¬ 
pute,  in  a  logical  manner,  the  bulk  modulus  of  the  sediment  1  water-mineral )  sy  stem  from 


its  components  without  estimations.  tins  subject  will  be  discussed  in  detail  in  t lie  next 
section  Prior  to  these  discussions,  the  pertinent  equations  of  elasticity  should  be 
reviewed.  The  equations  favored  are  those  involving  the  two  measured  constants  (den¬ 
sity  and  compressional-wave  velocity)  and  the  bulk  modulus. 

The  basic  equation  loi  the  velocity  of  a  compressional  wave,  I’  ,  is 


x  =  incompressibility  or  bulk  modulus  =  1/d 
p  =  shear  (rigidity)  modulus 
p  =  density 
i3  =  compi ’ssibility 

When  a  medium  lacks  rigidity,  p,  equation  1  becomes 


/  K  \  '  , .  1 

-  I  —  I  ,  or  K  =  p  r 

\p/  P 


V  =(— ) 

P  Id  p) 


Equation  2  applies  to  any  liquid,  emulsion,  or  suspension  which  lacks  rigidity. 

For  a  unit  volume  of  a  suspension  or  porous  material,  lacking  rigidity  and  com¬ 
posed  of  water  and  mineral  grains,  compressibility,  |J.  and  density,  p,  in  equation  2b  have 
been  expanded  into 


or,  expressed  as  k 


n  (k  -k  )  +  K 
'  s  w  w 


p  =  up. .  +  (l-n)p. 


n  -  volume  of  port  space  occupied  by  water  i .  iracnonai  porosity),  subscripts 
s "  and  “w"  indicate  mineral  solids  and  water,  jJK,  is  computed  with  equation  2a, 
and  (3S  with  the  Voigt-Reuss-Hill  averaging  method  i  Jiscussed  below );  and  "sic"  is 
used  to  indicate  system  moduli  computed  with  these  two  components  alone, 
and  to  differentiate  such  moduli  from  the  system  bulk  moduli,  k  fas 
in  equations  7- 10). 

The  result,  when  expanded  (3  and  p  (eqs.  3  and  4J  are  used  in  equation  2, 


p  ([-idK.  +  0")i3J][«PK.  +  <l_n)Pj]) 


is  known  as  the  Wood  equation  (Wood,  194 1 ),  its  experimental  and  theoretical  justit ica- 
tion  was  affirmed  by  Urick  ( 1947)  and  Chambre  ( 1955 ). 
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The  basic  equation  tor  the  velocity  ul  a  shear  wave,  1^,  is 

Substiuitmg  p  ~  p  ls~  into  equation  1  yields 


>h2 -*>■"?) 


which  is  the  best  equation  lor  the  detenu, n.ili.ui  ol  the  hulk  modulus,  x,  when  1,1, 

1  ft  s 

and  f)  are  known 

In  the  present  studs  .  i;1  and  p  were  measured,  and  x  computed  (as  discussed  in 
the  next  section)  These  three  constants  were  then  favored  in  computations  of  the  other 
constants  (when  was  not  known)  The  equations  used  were 


M  =  [pVr  "  K)  3  4 

-‘l  p 


Poisson's  ratio,  o  = 


-V;  -  pi'  2  (!'  I  I2  -  2 


Lame's  constant,  X  -  k-  2/3  p 


RESULTS  AND  CONCLUSIONS 
The  Bulk  Modulus  of  Saturated  Sediment 

INTRODUCTION 

The  best  values  tor  dynamic  bulk  moduli  ot  water-saturated  sediments  arc 
obtained  when  densities  arid  compressional-  and  shear-wave  velocities  are  know  n  (eq.  7). 
k  =  Up*  -  4 lj','3).  Lacking  shear-wave  velocities,  the  problem  is  to  compute  in  a 
logical  manner,  w  ithout  empirical  estimations,  or  constants,  values  for  bulk  moduli  which 
can  then  be  used  with  measured  densities  and  eompressional-wave  velocities  to  compute 
the  other  elastic  constants. 

As  noted  in  a  previous  section,  it  has  been  demonstrated  that  compressibilities 
and  their  reciprocals,  incompressibilities  or  bulk  moduli,  can  be  computed  for  a  unit 
volume  ol  a  rock  or  mineral  aggregate  Irom  the  volume  tric  contributions  cl  us  compo¬ 
nents.  The  question  is  can  an  aggregate  theory  be  appiied  to  relativelv  highly  porous 
saturated  sediments9  As  a  first  approach  to  this  question,  experimental  work  done  b\ 
Uriel;  1 1947  )  proved  that  an  aggregate  theory  in  the  form  of  equation  3a  could  be  applied 
to  suspensions  c.  ';aohnile  in  distilled  water.  Urick  used  a  deflocculent  in  these  experi¬ 
ments,  which  prev  ted  flocculation  of  the  kaolinitc  (  to  form  clay-mineral  structures), 
thus,  the  material  was  a  true  suspension  without  rigidity  and  the  Wood  equa¬ 
tion  (5 )  applied. 


DENSITY  (SOUND  V  E  IOC  I T  Y  |!  ,  dynei/t*’ 


Hamilton  ct  al.  ( computed  elastic  constants.  using  equations  Ml.  Ini 
natural  marine  sediments  using  measured  densities  and  compicssional-wave  velocities, 
and  values  tor  aggregate  compressibilities,  |J  .  computed  with  equation  3a.  except 
that  mineral  eoupressibilities  v.eie  computed  with  the  Reuss  method.  Such  computa¬ 
tions  imply  (see  cq.  1)  that  d  pip'  >  k,  all  the  excess  is  due  to  d>  uamie  ligidily  .  p, 
which  is  theoretical!),  present  (as  noted  by  Nate,  and  Drake,  l‘»i>3,  p.  SOS)  hut  apt  to 
be  vciy  small  Lacking  shear-wave  velocities  (<  ■  other  independently  measiiied  constants), 
this  method  could  not.  at  that  time,  be  evalu.i  cd 

T  here  are  now  sufficient  data  from  laboratory  and  field  experiments  wheie 
shear-wave  velocities  and  othei  constants  wcie  measured  to  allow  evaluation  ol  aggiegate 
theories,  as  applied  to  saturated  sediments,  and  to  improve  the  methods  lor  computing 
bulk  moduli.  In  the  cases  where  compressional-  and  shear-wave  velocity  and  density 
have  been  measured  in  a  salinated,  porous  sediment  or  rock,  a  true  value  of  the  sy  .tem 
bulk  modulus,  k  (eq  7).  is  always  greater  than  an  aggregate  bulk  modulus,  k  .  detei- 
mined  from  equations  3a,  b.  The  difference  lies  in  a  nccessaiy  additive  component  ol 
the  bulk  modulus  of  the  system:  lire  “skeletal,  structliial,  or  fiame  hulk  modulus." 

With  this  increment,  the  aggregate  theory  appeals  in  be  valid  tor  computations  of  bulk 
moduli  and  compressibilities  of  saturated  sediments. 

In  tigure  I .  p  f  •  is  plotted  against  an  aggregate  bulk  modulus  .  eq.  3n) 
which  does  not  include  the  frame-bulk  modulus.  If  the  Wood  equation  (5)  applied  tor 
these  deep-sea  sediments,  they  would  be  suspensions,  lacking  rigidity  and  frame-bulk 
moduli;  consequently,  the  Irne.pL  *  -  k,  shown  on  figure  I  would  adequately  define 
these  sediments.  The  tact  that  they  all  fall  above  the  line  indicates  that  k  is  not 
adequately  defined  or  that  rigidity  is  present,  or  both 


figure  I .  Hulk  modulus.  knw  i  without  a  flame  bulk  ir.snlulus,  -  ,,  1  ^  . 
jb\  seal  hill  (squares)  jn*i  jt>>  ssaJ  plain  I  triangles)  environments.- 


Il  has  become  apparent  (icstime  below)  lliat  almost  all  mantle  sediments 
have  measurable  ngnlitics  and  transmit  shear  waves  This  laet.  plus  the  need  lor  a 
name-bulk  modulus  as  an  increment  to  aggregate  hulk  moduli  (discussed  above), 
indicates  that  deviations  from  l he  Wood  equation  (as  m  tig.  1)  aie  due  to  the  presence 
ol  both  ngidtly  and  a  fume  hulls  modulus.  Pot  computations  ol  better  values  ol 
both  ngtdity  and  shear -wave  velocities  (and  otliei  constants),  the  frame-hulk  modulus 
has  io  he  considered  in  compulations  ol  the  hulk  modulus  ol  the  system. 


COMPONENTS  OF  THE  BULK  MODULUS 

lit  his  study  ol  the  elasticity  ol  porous  media  (rocks  and  sedimenls).  Uassmann 
I  I 5  1  i  considered  the  "closed  system"  in  which  the  water  within  pore  spaces  does  not 
"ciiculalc"  undei  the  small  stress  ol  an  elastic  wave,  and  the  elastic  moduli  ate  calcu¬ 
lated  with  the  equations  ol  elasticity.  Allhough  Biot  (  1941)  must  he  ciedilcd  with  the 
e.iilv  discussion  ol  the  "closed  sy  stem"  thorn  stal  ic  considerations),  (lassmaim  (  1 0  5  I ) 
has  a  particularly  clear  disenssion  ol  the  factors  involved  in  compulations  ol  dynamic 
elastic  constants  in  water-saturated,  porous  media 

Uassinann  (  |9S  I  )  demies1  rated  'll"  piw-ibibt;1  a :  ■  ’  uokjy.  ol  dci  ;>  lug  j  hulk 
modulus  tor  salinated  tucks  from  measniements  ol  eontpressional-  and  sheai-wave 
veloetlies  and  density  on  dry  material.  Ills  eqiulions  (discussed  below)  ate  used  m  a 
different  way,  in  this  report,  to  derive  an  essential  component  of  the  system  bulk 
modulus  (a  "frame  bulk  modulus"),  and  then  to  compute  bulk  moduli  lor  the  sedi¬ 
ments  i  i  this  report. 

In  dei tying  Ins  equation,  (.iassm.mn  used  the  principle  ot  el lect ive  pressure 
from  the  In'id  ol  soil  mechanics,  lot  details  see  any  leM  book  on  soil  mechanics 
(Skempton,  Idol,  has  a  particularly  good  discussion) 

)'  =  /’,  > /’  I  I  2) 

h  < 


where 


/’  =  (oral  piessuie 

/’ ,  =  hvdrostalic  pn  ssiire 
h  ' 

P  =  eflective.  or  intergranular,  pressure 

Upon  application  ol  a  total  pressure. />,  upon  all  sides  ol  a  unit  volume  ol  a 
porous,  "closed-system"  medium,  there  are  three  pressure  eltects.  (!)  the  ellecl  ol 
hydrostatic  picture,  /),  .  on  the  poie  water,  (2)  the  effect  of  hydrostatic  pressure,  P^  . 
on  the  mineral  grains  which  form  rhe  frame,  structure,  or  skeleton,  of  the  sediment, 
and  (3i  the  effective,  or  mtergranulai ,  ptesstire,/’  .  on  the  frame.  These  pressures, 
when  related  to  the  compressibilities  of  the  inineial  and  water  components  of  the 
unit  volume,  result  m  three  component  compressibilities  (ot  then  reciprocals,  bulk 
moduli).  At  an\  given  pressure,  these  bulk  moduli  are  dependent  on  porosity,  and 
are  related  to  the  system  bulk  modulus,  k,  by  tiassmann’s  basic  equations  tCiassmann, 
1951,  p.  15): 


s,+£> 


Q  = 


\-  +  C> 


II  ( K  -  K  ) 


(13a) 


wlic  re 


k  =  aggregate  bulk  modulus  of  mineral  solids 
K  —  bulk  modulus  ol  pme  watci 

k j  =  frame  bulk  modulus  ("skeletal''  bulk  modulus  of  Gassmanu) 
n  =  fractional  porosity. 

Gassmann’s  paper  (1951)  is  in  German  and  lie  did  not  show  the  complete 
derivation  of  equation  Ida;  therefore,  it  is  fully  developed  in  Appendix  A. 

When  solved  lory  equation  13a  becomes 


k/  n  (k  -  k.  )  +  K  I  ■  K  K 

K.  - - (  13b) 

J  n  (k  -  k  )  +  k  (k/k  ■  \) 

s  **  w  5 

Wh»n  the  bulk  modulus  of  the  skeleton  or  frame,  sj  ,  is  zero,  equation  13a 
reduces  to  the  Wood  equation  IS)  for  a  suspension  of  mineral  grains  in  water.  When 
porosity  is  zero  (r.e.,  a  solid  cube  of  minerals),  equations  13a  and  b  reduce  to 

K  =  K  =  K 

f  s 

White  and  Sengbush  (1953)  and  Wyllie,  Gregory,  and  Gardner  (1956),  with 
proper  reference  to  Gassmar.n  ( 1951),  came  to  the  same  conclusions  regarding  the 
relationships  of  the  bulk  moduli  as  shown  in  equation  13a.  Laughton's  (1957) 
“structural  bulk  modulus"  is,  ot  course,  fundamentally  the  same  as  Gassmann's 
skeletal  bulk  modulus. 

The  only  difference  between  Gassmann’s  skeletal  bulk  modulus,  k  ,  and 
Laughton’s  structural  bulk  modulus,  k  ,  is  that  in  the  fundamental  equations  relating 
pressures  and  compressibilities,  Gassmanu  included  effective  pressure  on  mineral 
compressibilities  (ie..  PfBs .  or  P{  ! k $),  and  Laughton  omitted  tins  effect  As  a  result, 
at  high  porosities,  k  ^  /  and  k  are  dose  together,  but  diverge  at  lower  porosities 
because,  as  porosity  approaches  zero,  hydrostatic  pressure  approaches  zero,  and  total 
pressure  approaches  effective  pressure  (eq.  12).  Laughton’s  equation  (x  =  Kf>(  +  K  , 
wherein  nm  is  k  in  equation  3b)  results  in  the  unacceptable  solution,  at  zero 
porosity,  that  k  =  k  ■  k£  Thus,  in  a  rigorous  solution,  Gassmann’s  prior  equations 
should  be  used,  especially  with  the  lower  porosities  of  sands  and  rocks. 

Using  Gassmann’s  equation  (13a)  to  derive  the  system  bulk  modulus  lor 
saturated  sediment  (without  a  gas  phase),  one  needs  to  know  four  components; 
porosity,  n,  the  bulk  modulus  of  the  pore  water,  k  ;  the  aggregate  bulk  modulus  ol  the 
mineral  solids,  k  \  and  the  frame  bulk  modulus,  k 

The  basic  assumption  underlying  computations  ol  the  bulk  modulus  ot  pore 
water  in  sediments  (in  this  report)  is  that  pore  water  and  bottom  water  have  the  same 
salinity.  Siever  cl  al.  (1965)  and  Friedman  cl  al.  (19p8)  have  shown  that  this  assump¬ 
tion  is  valid  within  small  limits.  Any  changes  ol  sound  velocity  and  density  caused 
by  these  salinitv  differences  can  be  disiegarded.  Thus,  the  values  ol  the  bulk  modulus 
of  distilled  and  seawater  can  be  found  in  tables  or  computed  with  equation  2a,  b  and 
used  as  the  bulk  modulus  of  sediment  pore  water,  .  Some  typical  values  for  sea¬ 
water  are  listed  in  Appendix  B. 
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In  iCl'CM  year*.  1  lie  elaMk  a»n*lants  ol  the  more  common  mmeiah  have  either 
been  Uciennincd  01  revised  h>  a  number  ol  specialists.  Aticnimn  is  especial l\  directed 
to  two  recent  compilations  which  also  include  excellent  discussion*  and  numerous  re  1c i - 
cnees  (Buell,  Anderson  and  i.iehcmiann.  l^OS).  because  ot  their  common 

x'cuirencc  in  marine  sediment'*,  altcntion  is  also  directed  to  reports  on  ealcite 
(I’csclmek.  I9t>2),  on  quart/  (McSkimn.  c/  a! .  ll)hS;  Soga,  ll)bh),  on  micas,  homhlcndc, 
and  leldspars  (Chnstensen.  ll,o5.  IVnn).  on  obsidian  (Manphnam  <7  ai,  I'Jtdi).  ami  the 
bulk  moduli  ot  oxides  (And;  i  .  ui i J  Nalc, 

The  availability  ot  these  ri\cn».  accurate  measurements  ot  the  elastic  constants 
ol  minerals  allows  use  ol  the  V  »npi  Reuss  l Irll  averaging  method  (heiealter  the  "VRH 
method"),  to  compute  the  bulk  modulus  id  a  rardomls -m lenlcd  mineral  aggregate  with 
some  eonliderice  (reference*  under  paragraph  nuinbcied  2.  above)  As  Hill  [  l*/52,  I^M) 
pointed  out,  the  Rcuss  average  assumes  uniform  stress  within  the  aggregate,  and  a  lower 
limit  to  ineoinpressibit its ,  k  is  computed  with 

A 


The  Voigt  averaging  method  assumes  uniform  strain  with  the  aggregate,  and  an  upper 
liurit  lu  incompressibility  ,  r.(  .  is  computed  with 


-'ik, 


where 


lj  and  I",  =  fractional  concentration  ol  phase,  (mineral  species)  in  a  unit 
volume 

K|  and  k  ,  =  incompressibilities  (hulk  moduli)  ot  the  phases 

Mill  showed  that  the  correct  average  lie,  hallway  between  the  Reus,  and  Voigt  aver¬ 
ages;  thus,  the  computed  bulk  modulus,  k,  ol  the  aggregate  is 

is  =  (  k  +  k  );  2 
A  I 

Mill  (l%3,  p.  3hl)  noted  that  both  the  Keuss  and  \’oigt  averages  weie  "rather  pom 
when  tbrt  nl' jse  nvjuuls  differ  bv  rvirc  thuu  u  fuctor  o!  2  or  vj."  T!'.is  i>  s!io\\u  !«*  be 
true,  lor  example,  if  the  VRH  method  is  tried  for  a  two-phase  suspension  composed  ot 
water  and  mineral  particles 

The  aggregate  bulk  modulus  of  mineral  grains,  s  ,  was  computed  in  this  study 
bv  using  the  VRH  averaging,  method. 

The  only  remaining  constant  needed  to  compute  the  bulk  modulus  of  a  satu¬ 
rated  sediment  (using  eq.  13a)  is  the  frame  bulk  modulus,  .  The  recommended  way 
to  dense  this  value  (for  use  in  eq  13a)  is  to  enter  a  curve  feinting  this  constant  to 
porosity  or  void  ratio  (volume  ol  voids/volume  ol  solids).  These  curves  have  been 
derived  in  two  ways  by  Laughton  (  1957)  using  drained,  static  complexion  tests,  and 
in  tins  report  by  solving  equation  ( 13b)  with  carelully  selected  datj  which  include 
values  lot  density,  and  compressional-  and  sheat-wave  velocities,  as  described  below. 

Laughton  t  1957)  derived  two  curves  lelatmg  static  ' r.itne  bulk  moduli  to  void 
ratios  through  a  series  of  drained  compression  tests  on  two  deep-sea  sediment  types 
a  calcareous  ooze  and  a  “terrigenous  mud  ”  As  Laughton  noted,  a  static  tc !.  because 
of  various  laboratorv  testing  difficulties,  has  wide  margins  ot  error  (he  estimated  plus 
or  minus  2,5  X  10ll>  d\ nes/cm' for  Ins  tes's).  Although  theoretically  usable  lor 
dynamic  compulations,  a  static  compression  lest  loi  the  I  tame  bulk  modulus  is  less  upl 


l'1  be  as  collect  as  a  dynamically  dclcnnincd  modulus,  cspcviaih  loi  (lie  puipoccs  ol 
computing  olhci  dynamic  claslic  constant  winch  aic  usaldc  in  tiiidciwalci  sound  oi 
goopliy  sical  studies. 

DYNAMIC  I- KAMI:  HULK  MODULUS 

A  "dynamic  liame  bulk  niinlulus."  Sj  .  as  used  in  tins  icpoit,  is  one  detei- 
mined  with  equation  l.'b.  using  dynamic  values  lei  (he  saiious  bulk  moduli  icquucd  to 
solve  the  equation. 

1  be  accumulation  ol  nicasuicinculs  ol  the  elastic  constants  ol  mmctals.  and 
compiossioiiai-  and  shcai-wavc  velocities  and  densities  in  locks  and  sediments  during 
the  past  15  yeais.  now  permits,  in  the  author's  opinion,  computations  id  s\  stent-hulk 
moduli  (k  in  cqs.  7  and  Ida),  and  values  loi  the  hulk  moduli  ol  poie  watet  and 
m, metals  (x  and  x5 1  loi  enough  locks  and  sedimcnls  to  dclinc,  approximately,  [elation- 
slops  between  dynamic  frame-bulk  moduli  and  poiostty  ( s(  ami  n  in  eq.  Ida,  h). 

Data  winch  can  he  used  tit  computations  ol  a  dynamic  Itanie-lmlk  modulus, 
using  equation  Idb,  imisl  meet  the  hallowing  requirements. 

1.  Measurements  must  include  both  Cc'iuptessional-  and  slieaiwjve  velocities  in 
walei-satiiiuled  sediments  oi  nicks. 

2.  Density  and  porosity  are  known,  oi  can  he  computed  from  known  data. 

d.  The  r  lineralogy  of  the  medium  is  known  (including  peicentage  of  volume 
oi  species),  oi  can  he  reasonable  estimated. 

4.  The  nature  of  the  poie  fluid  is  known,  oi  can  he  reasonably  mteiied. 

5.  The  lempeiatutc  within  the  medium  at  the  tune  ol  the  velocity  ineasuic- 
nicnts  is  known,  or  can  be  reasonable  inferred. 

f>.  The  pressure  on  the  sample,  or  medium,  ;s  known  or  can  he  computed 
(because  any  bulk  modulus  varies  with  pressure). 

A  selection  of  data  (Ironi  (lie  literature  and  Irom  field  work  by  the  author 
and  Bucket)  which  met  the  above  requirements  was  used  to  compute  dynamic  flame- 
hulk  moduli  for  a  variety  ot  saturated  rocks  and  sediments,  these  moduli  were  then 
plotted  against  porosity  (fig.  7).  Ail  data  in  figure  2  wore  referted  to  an  approximate, 
common  pressure  which  would  be  expectable  with  1  motet  oi  less  ol  oveibuidcu  pres¬ 
sure,  and  a  common  temperature  of  2d°C;  tempeiaturc  is  particulars  important  in 
computing  compressibilities  of  pore  fluids. 

In  a  previous  section  it  was  noted  that  in  Gasstnann's  equation  (l.'b),  when 
poiusilv  is  zero,  the  bulk  modulus  ot  the  frame  equals  the  bulk  modulus  ol  the  solid 
mineral  (or  the  aggregate  of  minerals),  this  allows  a  valid  point  to  be  plotted  at  zero 
porosity  through  whitli  any  regression  line  lot  media  of  the  same  mineralogy  should 
pass.  Four  such  points  are  noted  at  zero  porosity  in  figure  2:  quartz  (Soga,  1968). 
calcile  (I’eselniek,  1962),  and  two  aggregate  values  for  sediments  off  San  Diego 
(Appendix  B) 

A  single  regression  line  is  shown  in  figure  2  for  all  dala.  However,  this  line 
should  be  considered  as  approximate  because  of  the  sparse  data  and  because  there  should 
be  individual  lines  for  each  sediment-rock  nuneial  type,  Lxanipies  oi  these  lines  should 
he  I  I )  relatively  pure  calcine  limestones,  sands,  anu  calcareous  oozes,  (2)  relatively  pure 


FRAME  BULK  MODULUS,  K, .  *  10  dynej/cm 


1000  0 


I\» 


CALClTE 


1 - 1 - T 


QUARTZ 

\  * 


100  0 


10.0 


1  0 


0  1 


0.01 


•  v  • 


i — r 


\s. 


■  LIMESTONE 
•  SANDSTONE 
O  SAND,  IN  SITU 
Ol  SAND,  LAB 
A  SHALE 

A  CLAY. SILT,  IN  SITU 
ix  STATIC  (LAUGHTON,  1957' 


1 - 1 - 1  I 


\ 


\ 


A. 


\ 


10 


20 


30 


40  50  60 

POROSITY,  % 


70 


80 


90 


100 


I  rgure  ,\  l\uo*at>  u;.n:.n.i,  iiamv  L..IL  ranjuko,  „ ,.  X:c  test  for  liisoiwiuii, 
vc  Append  is  H  ! r  regies, ion  -.-qua  I  ions  Let  .,i iul >  jiul  sill -etas  s. 


quail/  sandstones  and  sands,  and  (5)  naniijl,  null  I  imnioiaiic  scdimentaiy  rocks,  sands, 
and  finer-grained  teriigemju.s  sediments.  However,  il  is  surprising  how  well  tiro 
generalized  ngression  line  defines  the  available  data,  for  example,  lire  trend  ol  the  dau 
liont  7cru  porosity  to  ahoul  45  percent  poiosiiy  would  have  predicted  reasonable  values 
ol  x  f  for  the  high  porosities  above  70  percent  Appendix  L  provides  regression  equa¬ 
tions  lor  sands  and  silt-class. 

Because  the  frame-bulk  modulus  vaiies  with  porosity  and  effective  pressure 
(Appendix  A),  a  line  relating  log  Kf  and  porosity  tin  pure  caleitiv  locks  and  sediments 
should  he  above  those  f.  i  pure  quartz  and  normal,  inulliinineral  earth  materials,  as 
porosity  incieases  from  zeu>.  This  is  indicated  i it  figure  2  by  the  values  lor  ealeite  at 
zero  porosity,  limestone  (l’esehnck,  1962).  and  Austin  Clijlk  (White  and  Senghush, 
1953).  eoinpaied  with  the  values  lor  quartz  (Soga.  196, H)  and  quaitz  sands  (Shell 
Development  (’>>.,  private  enmm.,  1965)  In  Laughton’s  static  tests,  the  line  relating 
Kj  with  vo'i  ratio  for  calcareous  ocz”  was  above  and  roughly  parallel  to  that  for 
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lemgcnous  mud  (195/.,  I'm.  14.  p.  25$).  The  lowest  pressures  at  which  Laughton 
obtained  static  measurements  of  frame-bulk  moduli  (Ins  x  )  calcareous  ooze  (an  esti¬ 
mate  at  P  =  S  kg  em: )  and  terrigenous  mud  (P  =  lo  kg.  cm^  )  are  included  in  figure  2 
lor  comparison  with  the  dynamic  moduli.  Other  sources  of  data  for  points  in  figure  2 
aie  sandstone  (king,  19(->(>),  Lagle  l  ord  shale  (White  and  Sengbnsh.  1953).  and  sands 
and  silts  clays  oil  San  Diego  (Hamilton  cl  j I9ti9). 

Computations  and  Discussions  of  Elastic 
Constants 

INTRODUCTION 

The  elastic  constants  lor  various  sediment  ty  pes  w  ithin  the  three  environments 
were  computed  with  equations  7  to  II.  The  input  data  were  (  I )  measured  values  of 
saturated  bulk  densities  and  compressional-wave  velocities  (Part  I,  tables  I  and  2:  re¬ 
produced  in  this  report  as  tables  0-1  and  D-2,  Appendix  D)  and  (2):  computed  values 
of  bulk  moduli,  k,  using  equation  Ida  with:  (at  measured  values  of  porosities,  n 
(Part  I.  tables  1  and  2),  (b)  computed  values  lor  the  bulk  moduli  of  the  pore  water, 

,  using  equation  2a,b.  (c)  computed  values  for  the  aggregate  bull,  moduli  of  the 
mineral  grains,  ks,  using  literature  values  for  ihe  bulk  moduli  of  individual  minerals 
and  summing  their  individual  volumetric  contributions  to  the  aggregate  bulk  moduli 
using  the  VRI1  method:  and  (d|  values  computed  for  the  frame-bulk  moduli,  k  .  by 
entering  the  regression  equations  for  sands  and  clay -silt  (fie.  2,  Appendix  ('). 

The  computed  elastic  constants  are  listed  in  tables  I  and  2.  The  relation¬ 
ships  of  several  constants  with  other  physical  properties  are  shown  in  the  figures. 

All  of  ihe  mput  data  noted  above  have  margins  ot  error,  some  known  and 
some  unknown.  All  of  the  computed  constants  should  be  considered  approximations. 
Consequently,  no  attempt  was  made  to  statistically  estimate  variances  or  errors  in  the 
final  computations.  The  numbers  of  decimal  places  shown  in  the  tables  are  for  purposes 
of  comparison  between  Ihe  various  computations  ar,d  should  not  be  taken  as  the 
author's  estimate  of  accuracy. 

Computations  of  elastic  constants  are  strongly  dependent  on  accurately 
"■easurei  values  of  density,  puiosity.  and  compressionai  velocity.  An  examination  of 
equation;  I  and  Ida  indicate  why  this  is  true.  Ifpl  -  is  not  greater  than  the  com¬ 
puted  bulk  modulus,  k  (which  is  strongly  dependent'  on  perosity).  there  is  no  rigidity 
or  shear  wave,  and  Pott-on ‘s  ratio  is  0.50.  Because  the  rigidity  modulus,  p,  is  so  small 
in  most  marine  sediments,  inaccurate  values  of  p.  L  .  and  n.  especially,  may  resuit  in 
computations  indicating  no,  or  too  great,  rigidity.  £or  these  reasons,  all  attempts 
should  be  made  in  the  laboratory  to  improve  measurements  of  these  properties.  For 
example,  "saturated  densities"  should  be  truly,  fully  saturated  (not  always  easy  in 
sands),  and  properties  should  be  "salt-free"  j i.e:,  no  dried  salt  weighed  with  the  dried 
minerals,  see  Part  1).  Accurate  measurement  of  bulk  gram  density,  without  salt,  aids  in 
cross-checking  density  and  porosity  values  (Part  I).  In  the  computed  elastic  constants 
listed  in  tables  1  und  2,  k  >  p'lp~  m  some  principal  sediment  types,  as  follows  3  out  of 
of  40  cases  in  continental-terrace  clayey  silts.  1  out  of  52  cases  in  abyssal-plain  silt- 
clavs,  none  out  of  41  in  aby  ssal-hill  silt-clays,  and  2  out  of  9  in  fine  sand,  these  cases 
were  not  included  in  the  averaces  of  p,  i  .  ot  o.  For  this  reason,  but  mostly  because  ot 

j  • 

rounding  oti  at  a  lower  number  of  decimal  places  than  used  in  the  computations,  the 
average  values  in  the  tables  cannot  usually  be  exactly  interrelated  with  the  elastic 
equations. 

It  should  be  emphasized  that  all  values  for  the  elastic  constants  aie  for  a  temp¬ 
erature  of  23°C  and  I  atmosphere  ot  pressure.  Methods  for  relating  these  constants  to 
in  situ  values  will  be  discussed  in  Pari  III  (Prediction). 


TABLE  I :  CONTINENTAL  TERRACE  (SHELF  AND  SLOPE)  ENVIRONMENT:  COMPUTED 
ELASTIC  CONSTANTS  IN  SEDIMENTS 


I 

I 

i 

J 

) 


I.ame's  constant 


EFFECTS  OF  USING  A  FRAME-BULK  MODULUS 


A  value  tor  a  “system  bulk  modulus,”  x,  can  be  computed  from  its  compo¬ 
nents  (and  porosity)  in  two  ways.  ( 1 )  using  the  moduli  for  pore  water,  xw.,  and 
minerals,  x5.  in  equation  3b,  and  assuming  that  xJlt.  is  x,  or  (2)  using  moduli  for  pore 
water,  minerals,  and  a  frame-bulk  modulus,  Kj,  in  equation  13a.  The  effect  of  using  a 
fiame-bulk  modulus  is  to  increase  the  value  of  the  system-bulk  modulus  by  a  small,  but 
significant  amount.  When  measured  values  of  density  and  comprcssional-wave  velocity- 
are  used  with  computed  values  for  x  to  compute  rigidity  ,p  (eq.  8),  an  mciease  in  x 
reduces  the  value  of  rigidity  .  Lower  rigidities  cause  lower  shear-wave  velocities  (eq.  6). 
Numerical  examples  of  the  ufferences  in  x  ,  p,  and  I  5  when  x  is  computed  with  and 
without  x.  are  given  below. 


System  Bulk  Modulus  Used* 

( dynes  eni-  X  1 01  u  ) 

U 

(dynes. cm-  X  1010)  ** 

>'r 

(msec)*** 

1.  Fine  sand  (>i=4 Srr) 

x JK,( without  K,.)  5.0333 

0.6997 

594 

k  (with  x y)'.  5.5418 

0.3183 

401 

2.  Abyssal-hall  silty  day  (fl=84.3‘e> 

ksw  (without  Kj).  2.8165 

0.1151 

296 

x  (with  Kf)\  2.8292 

0.1056 

284 

*K„  eq  5b;  K.  eq.  13a 
svv  ^ 

“M  =  Npf  p2  -  K)3,4 
•**f  s  =  ij \jjp)xa 

Differences  in  rigidities  and  shear-wave  velocities  are  greater  in  the  lower- 
poiosity  sands  (above)  because  x^  increases  with  decreasing  porosity  (fig.  2),  causing 
a  larger  difference  between  the  two  bulk  moduli  (xSH  and  x). 

COMPRESSIBILITY  AND  INCOMPRESSIBILITY  (BULK  MODULUS) 

The  theoretical  and  experimental  evidence  concerning  compressibility  and 
the  bulk  modulus  have  been  discussed  in  previous  sections;  therefore,  this  section  will 
be  confined  to  empirical  relationships. 

Because  the  bulk  modulus  is  favored  in  the  computations,  values  of  compressi¬ 
bility  are  not  listed  in  tables  1  and  2.  However,  if  desired,  they  can  be  easily  computed 
ld='l/Kl. 

At  am  given  temperature  and  pressure,  values  of  compressibility  and  the  bulk 
modulus  vary-  within  a  small  range  in  both  pore  water  and  mineral  solids  (Appendix  B). 
The  volumetric  contributions  of  both  of  these  moduli,  and  the  frame-bulk  modulus,  are 
strongly  dependent  on  porosity,  n ,  space  occupied  by  water,  or  (1  •  .n),  the  space 
occupied  by  mineral  solids  in  a  unit  volume.  Therefoie,  it  is  no  surprise  to  find  that  the 
plots  of  computed  values  of  j3  or  x  vs.  n,  or  density,)},  show  a  good  correlation 
(figs.  3-10). 
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f  igure  5.  Density  vs.  compressibility . 
continental  terrace. 
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BULK  MODULUS.  dynti/tn’  "  10'“  COMPRESSIBILITY.  dyn«i/cm’ 
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1’igurc  6.  Density  w  compressibility,  abyssal 
hill  (squares!  and  abyssal  plain  (triangles) 
environments. 
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f  igure  8.  Porosity  vs.  bulk  modulus, 
abyssal  hill  (squares)  and  abyssal  plain 
(triangles)  environments. 


6  80 
660 
6.40* 

j 

6  20  j 
6  00 1 
5  801 
5  60  j- 
5  40r 
5  20, 
5  00 1 


>'>  4,80; 


4  60 
4  40 
4  20 
4  OOj 
3  80| 


3  40 
3  20 

3  OOf 
20 


* 

\ 


V 

\ 


30  40  50  60  70  80  «0 
POROSITY,  % 


J 

100 
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BULK  MODULUS,  dynes./cm'  x  JO" 
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I-igure  V.  Density  vs.  bulk  modulus 
continental  lettuce. 
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I  igure  10-  Density  vs.  bulk  modulus,  abyssal  hill 
t square-.)  and  abyssal  plain  t triangles y  environments. 


Values  lor  2  o[  k  at  23°C  and  I  atmosphere  can  be  derived  by  entering  figures  3 
to  10  with  density  or  porosity,  or  by  use  of  the  regression  equations  for  the  bulk  modulus 
shown  in  the  figures  (see  Appendix  C). 

The  bulk  moduli  of  the  deep-water  sediments  are  plotted  against  pfp2  in  fig¬ 
ure  1 1 ,  together  with  a  line,  p  1  =  k  .  in  this  figure  the  values  of  «  contain  the  frame- 

bulk  modulus,  and  are  assumed  to  represml  “true”  values  of  k  (as  compared  with  fig.  I , 
in  which  kjk,  lacked  frame-bulk  moduli).  If  these  are  true  values  of  k,  the  divergence  of 
the  data  points  from  the  line  indicate  the  presence  (and  values)  of  rigidity  (fil'p  r  k  + 
4/3p). 
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RIGIDITS  (SHEAR)  MODULUS  AND  SHEAR  WAVE  VELOCITY 
INTRODUCTION 


I  here  lias  long  been  a  question  in  marine  gcopln  sics  and  underwater  acoustics: 
do  surficial,  natural,  marine  sediments  have  enough  rigidity  to  allow  transmission  of 
shear  waves;  and,  if  so,  what  arc  reasonable  values  of  these  moduli'’  The  answer 
appears  to  be  that  ( I)  in  some  bays  and  estuaries,  near  river  deltas,  and  in  a  few  other 
localities,  sediments  are  deposited  at  a  fast  rale,  lack  appreciable  structural  strength, 
and  for  all  practical  purposes  are  little  nioic  than  suspensions,  however,  (2)  almost  all 
of  the  remainder  of  continental-terrace  and  deep-sea  sediments  should  possess  enough 
rigidity  to  allow  transmission  of  shear  waves. 

The  evidence  for  the  above  statements  can  be  divided  into  experimental  evi 
dcnce  in  the  laboratory  and  in  situ  on  land  and  on  the  sea  floor,  and  theoretical  evi¬ 
dence.  Because  the  presence  of  a  shear  wave  is  dependent  on  rigidity  (eq.  6),  and  most 
measurements  have  been  of  shear  waves,  lire  following  discussion  wiii  be  largely  confined 
to  shear  waves.  Values  of  shear-wave  velocities  from  the  discussions  below  arc  assembled 
in  tabic  3. 


LABORATORY  MEASUREMENTS  OF  SHEAR  WAVE  VELOCITIES 


In  the  laboratory,  shear  waves  can  be  transmitted  through  dry  and  saturated 
sands  at  relatively  low  normal  pressures.  Experimental  difficulties  have,  so  far,  prevented 
low-pressure,  laboratory  measurements  of  shear-wave  velocities  in  high -porosity,  natural 
marine  sediments.  For  example,  Laughton  ( 1 957)  was  unable  to  measure  shear  waves  in 
a  sample  of  calcareous  ooze  until  the  material  was  compacted  to  a  porosity  of  about 
33  percent  under  a  pressure  of  5 12  kg, cm*. 

Much  of  the  laboratory  work  on  shear  waves  in  sands  has  come  from  soil- 
mechanics  and  oil -industry  research.  Recent  papers  which  summarize  (and  add  to)  much 
of  this  work  arc  by  Barkan  ( 1 962),  Hardin  and  Richart  ( 1 963,  including  discussions), 
and  Shell  Development  Company  (private  communication,  1965,  figs.  12  and  13,  this 
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TABU  3.  SUMMARY  01-  SHLAR-WAVL  VLI.OCITILS  IN 
FULLY-SATURATLL)  SFDIMKNTS. 


Sediment 

>* 

(m,  sec) 

Ref* 

Remarks 

Laboratory 

Coaisc  sand 

95 

1 

P-0.01  kg/cm ^ ;  resonant  col.;  Ottawa  sand 

Coaise  sand 

133 

1 

P  -  0.1 2  kg/cm2;  resonant  col.;  crushed  sand 

Coarse  sand 

285 

2 

P  =  1  4  kg  cm* ,  pulse  tech.;  St.  Peters  sand 

Fine  sand 

260 

2 

P  =  1 .4  kg/env;  pulse  tech.;  St.  Peters  sand 

ir,  Situ  (Land) 

Sand  (dense) 

250 

3 

Medium  sand 

160 

3 

•  Russian  soil-mechanics  lit. 

Fine  sand 

110 

3 

Clay 

150 

3 

• 

Sand 

130 

3 

Russian  soil-mechanics  lit.;  vibration:  10-35  Hz. 

Sand 

534 

4 

50-ft  depth;  explosive 

Clay-silt 

244 

4 

10-ft  depth;  explosive 

Clay-silt 

131 

5 

Sediment  sfc,  vibration  35  400  Hz; 

Fine  sand 

190 

6 

Sand 

255 

6 

Sand 

230 

6 

Sand 

227 

6 

Sediment  sfc;  vibration:  8-60  Hz 

Clay-silt 

150 

6 

...-silt 

180 

6 

Clay-silt 

200 

6 

Clay 

230 

7 

Bore  hole;  3-m  deptn;  explosive  (detonator) 

Clay 

120 

8 

0  to  20  m  depth  (seismic  logging);  from  i  u/l/s 

Clay 

170 

8 

0  to  10  m  depth  (seismic  logging);  from 

TABl.L  3  (continued ) 


Sediment  j 

(in/sec) 

Ref.* 

Remarks 

In  Sim  (Sea  Moot  I.  computed  from 

Stoncles  waves 

Sil  1  <1  ay 

50-190 

9 

Deep  Indian  Ocean,  explosive 

Sil  l-cla\ 

101 

1  i 

Sill-clay 

107 

10 

^  San  Diego  Trough,  explosive 

Sill-clay 

87 

10 

Medium  sand 

197 

1 1 

Com.  shell  off  San  Diego,  explosive 

Fine  sand 

101 

1 1 

Theoretical 

Pacific  Basin  j 

250 

12 

Pacific  Basin 

250 

13 

Atlantic  Basin 

150 

14 

.  Rayleigh  wave  model  (in  situ  values) 

Atlantic  Basin 

30-120 

15 

Argentine 

Basin 

200-900 

16 

j 

Computed  average  values,  this  report  (23°C;  1  atmosphere) 


Fine  sand 

382 

Continental  shelf 

Silty  clay 

28 1 

Continental  shelf 

Silty  clay 

232 

Abyssal  hills,  Pacific  Basin 

Clay 

195 

Abyssal  hills.  Pacific  Basin 

Silty  clay 

240 

Deep  abyssal  plains,  adjacent  to  Pacific  Basin 

•Refeio'cs 

1.  Hardin  and  Richart  (1963;  and  discussions) 

2.  Shell  Development  Co.,  personal  communication,  1965) 

3.  Barkan  (196?.) 

4.  White  and  Sengfcush  ( 1953) 

5.  Jones  (1958) 

6.  Heukelom  (1961) 

7.  Molotova  (1  96  3) 

8.  Zhadin,  V.  V.  {in  Molotova  and  Vassil’ev,  1960) 


9.  Davies  (1965) 

10.  Hamilton  a  al.  ( 1969) 

11.  Bucket  el  ai  (1964) 

1 2.  Kovach  and  Press  (1961) 

13.  Oliver  and  Dorman  <1961) 

14.  Latham  and  Sutton  <  196t>) 

15.  Andvisunand  Latham  ( 1969) 

16.  S>  kes  and  Oliver  ( 1964) 


report).  Many  ul  the  studies  referenced  rn  these  reports  also  include  values  lor  the 
velocity  ol  eompressional  waves  m  artificial  sands  at  “full  saturation''  (usually  with  dis¬ 
tilled  water)  Most  of  these  values  must  be  disregarded.  The  soils  engineer  must  know  the 
physical  properties  of  soils  from  the  dry  to  fully-satuiated  state,  and  the  author  does  not 
question  those  values  at  less  than  lull  saturation.  1  lowevcr.  unless  extreme  care  is  taken 
to  evacuate  all  air  from  the  “saturated"  samples,  the  values  for  eompressional -wave 
velocities  arc  tar  too  low.  Any  values  toi  eompressional -w  ave  velocities  at  1 00  percent  sale 
ration  should  be  well  above  velocities  in  water.  Sc h on  ( 1963)  illustrated  the  marked 
increase  in  eontprcssional-wavc  velocity  as  saturation  nears  UK)  percent.  Examples  of 
reasonable  values  lot  eompressional  velocities  in  artificial  sands  at  lull  saturation  are  in 
reports  by  tli.mdl  ( I960),  Schorl  (1963),  Nolle  cl  at.  (196.3),  and  Shell  Developm.’nt  Co. 
(see  tigs  13  and  13).  These  values  to;  artificial  sands  arc  consistent  with  those  measured 
in  natural  sands  m  the  laboratory  ,  and  in  situ  by  Hamilton  cl  ai  ( 1 956),  Sliumwav  (  1 960), 
and  with  the  measurements  ot  this  icport.  Tire  velocity  ol  sircar  waves  in  fulls  -saturated 
sands  is  slightly  lower  than  in  dry  oi  partially  saturated  sands  (Hardin  and  Kichart,  1963, 
Shell  Development  Co.,  figs.  13  and  13). 

Wilson  and  his  students  at  the  U.  S.  Naval  Postgraduate  School  have  a  continuing 
program  to  measure  dynamic  rigidity  in  artificially  sedimented  clays  and  natural  marine 
sediments,  using  a  torsional  oscillator  resonance  technique  (Hutchins  and  W  ilson.  1968. 
Hutchins,  1967.  Cohen,  1968).  Ty  pical  values  of  rigidity  have  been.  2.16  dynes/cni*  X 
10'  tor  a  32-perccnt  concentration  (by  weight)  of  kaolimtc  in  distilled  water,  and 
4.6  dynes/cm-  X  10^  for  a  1 9-pcrccn t  concentration  of  bentonite  in  distilled  w  ater 
(Cohen,  I9o8).  Using  equation  6  and  values  for  density  and  rigidity  troin  Cohen’s  report, 
the  computed  values  of  shear-wave  velocities  in  four  selected  samples  were  5  to  7  m/sec. 
Some  additional  aspects  of  (his  work  will  be  discussed  below. 

IS  SITU  MEASUREMENTS  OF  SHEAR-WAVE  VELOCITIES  (LAND) 

In  situ  experimental  wo'k  on  shear-wave  velocities  on  land  has  largely  been  by 
oil-industry  and  soil-mechanics  laboratories.  White  and  Scngbush  ( 1953)  measured  sircar 
and  eompressional  velocities  in  saturated  sands  and  other  materials  using  explosive  charges 
(blasting  caps]  and  dropped  weights  in  a  borehole.  Heukelom  (1961 )  used  dy  itarnic  deflec¬ 
tion  techniques  to  measure  rigidity  and  shear-wave  velocities  in  a  variety  of  natural  sedi¬ 
ments  in  the  Netherlands.  Russian  work  in  soil  mechanics  in  this  field  was  assembled  by 
Barkan  (1962).  Jones  ( 1958)  has  described  measurements  of  dynamic  properties  of  soils 
at  the  RoaC  Research  Laboratory  ,  Harmond.  worth,  Eng'and,  in  which  an  electromagnetic 
vibrator  was  used  U)  generate  elastic  waves.  Thcs:  measurements  are  of  particular  inter¬ 
est  because  tile  technique  and  theory  were  similar  lo  those  of  Davies  ( 1 9o5)  and  of  Bucket 
cl  ai  ( 1 964)  tor  in  situ  measurements  on  the  sea  floor. 

IS  SITU  MEASUREMENTS  OF  SHEAR  W  AVE  VELOCITIES  (SEA  FLOOR) 

Davies  (1965),  Bucker  cl  at.  (1964),  and  Bucket  and  Ins  colleagues  (in  Hamilton 
et  ai,  1969)  have  measured  values  of  Stoncley  -wave  velocities  front  which  shear  wave 
velocities  can  be  computed.  Davies’  measurements  were  in  the  deep  Indian  Ocean,  and 
Bucher's  were  in  sands  and  clay -sills  of  the  continental  shell.  The  conversion  ol  com¬ 
pressions)  waves  to  shear  waves  at  refraction  boundaries  within  decp-sca  sediments  was 
reported  by  Nate  and  Drake  ( 1957). 


THEORETICAL  VALU1  S  OF  SHEAR-WAVE  VtLOCITIES 


The  theoretical  values  of  shear-wave  velocities  noted  in  this  section  come  from 
the  reconciliation  of  experimental  data  with  theory  in  connection  with  the  propagation 
of  Love  and  Rayleigh  waves  across  ocean  basins,  and  from  the  computations  of  this 
report. 

To  get  quantitative  agreement  between  theory  and  observations  of  the  dispersion 
of  Love  and  Rayleigh  waves  across  ocean  basins,  Kov,  ch  and  Press  (1961),  and  Oliver  and 
Dorman  (1961),  used  shear-wave  velocities  of  250  m/sec  in  their  models.  Latham  and 
Sutton  ( 1966),  for  an  area  near  Bermuda,  used  a  value  of  1 50  m/sec,  and  Sykes  and 
Oliver  (1964)  used  200  to  400  m/sec  for  the  Argentine  Basin. 

Anderson  and  Latham  (1969)  studied  the  dispersion  of  Rayleigh  waves  caused  by 
sediment  layers  in  the  sea  floor  between  the  Mid- Atlantic  Ridge  and  Bermuda.  They 
derived  values  of  shear-wave  velocities  for  various  cases.  Because  sediment  properties  m 
the  area  are  reasonably  well  known,  a  value  for  shear-wave  velocity  in  the  upper  1  meter 
of  sediments  was  computed  using  the  method  described  in  this  report.  All  values  were 
corrected  to  in  situ  values  using  the  methods  of  Pari  111.  The  necessary  data  for  the 
sediments  were  taken  from  Horn  el  al.  ( 1 968c>:  average,  laboratory  values  of  density, 
porosity  ,  and  velocity  for  the  upper  1  meter  of  five  cores  taken  on  the  Mid-Atlantic  Ridge 
and  Sohm  Abyssal  Plain  (Stations  AS7-1 ,2,3,4,  ano  5)  weie  averaged  and  corrected  to 
in  situ  values.  The  velocity  of  bottom  water  at  5000  m  ( 1  54S  m/sec)  was  taken  from 
NAVOCEANO  TR-17'  ( 1965);  bottom-water  density  was  assumed  to  be  1.1)5  g/cc.  The 
computed  in  situ  value  for  the  average  shear-wave  velocity  ( 148  m/sec)  lies  between 
Anderson  and  Latham  Cases  C'  and  D  ( 160  and  70  m/sec).  The  authors  prefer  Case  D,  but 
cannot  climinaie  Case  F  which  corresponds  most  nearly  to  seismic  profiler  results  and 
calls  for  shear-wave  velocities  of  30  m/sec  at  the  lop  and  1 20  m/sec  at  the  bottom  of  a 
layer  1  50  n.  thick.  Similar  computations  for  a  generalised  red  clay  in  the  Pacific  Basin 
at  a  water  depth  of  5000  m  indicate  a  jcdiment-surface  shear-wave  velocity  of  about  220 
m/sec,  which  can  be  compared  with  rhe  value  of  250  m/sec  used  by  Kovach  and  Press 
(1961 1  and  Oliver  and  Dorman  ( 1961)  to  reconcile  experiment  with  theory  in  their 
studies  of  Rayieigh-wave  dispersion. 

Computed  values  of  shear-wave  velocities  for  each  sediment  type  within  each 
environment  in  this  report  are  listed  in  tables  I  and  2.  For  comparison  with  the  values 
from  the  studies  noted  above,  average  laboratory  values  for  several  Major  sediment  types 
are  also  sh  wn  in  table  3.  The  cotrection  of  laboratory  values  of  sediment  mass  physical 
properties  to  in  situ  values  is  a  subject  in  Part  III  (TP  145).  In  the  case  of  shear-wave 
velocity  in  surficial  sediments,  such  computations  indicate  a  progressive  decrease  in 
velocity  with  water  depth.  This  velocity  deurease  from  laboratory  to  5000  m  seawater 
depth,  for  a  high-porosity  silty  clay,  is  of  the  order  of  10  to  15  percent. 

An  examination  of  tabic  3  indicates  that  values  computed  in  this  report  are 
withun  the  range  of  (and  consistent  with)  a  wide  range  of  laboratory  and  in  situ  measure¬ 
ments  of  shear  waves  in  saturated  sediments. 

CAUSES  Ol  RIGIDITY 

In  saturated  sediments,  rigidity  is  related  to  sediment  structure  and  the  complex 
factors  restricting  relative  inlerparticle  movements  under  shear  stresses.  In  Part  I,  the 
common  sediment  structures  (Part  I.  fig.  1)  were  reviewed;  in  this  section,  those  aspects 
uf  sediment  structural  strength  relating  to  rigidity  will  be  discussed. 

Shear  strength  is  one  of  the  critical  engineering  properties  of  sediments  or  soils. 
As  a  result,  there  is  a  voluminous  literature  on  this  subject  in  the  fields  of  soil  mechanics 
and  foundation  engineering.  Fully  refeier.ced  discussions  appear  in  recent  textbooks 
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( e.g Jumikis,  1964:  vong  and  Warkentin,  1966)  and  in  recent  papers  of  particular 
interest  to  'he  present  discussion  (Rosenqvist,  1960;  Lambe,  1960;  Schmertniann  and 
Osterberg,  I960;  Hv  aislev,  1960;  Seed  and  Chan,  1961 ;  Mitchell,  l964;Grim,  1962). 
The  following  outline  is  derived  from  these  sources  and  those  referenced  below. 

The  shear  strength  of  a  sediment  at  failure,  t.-.  is  represented  by 

Tf  =  r  +  (o  -  p)  tan  c?  (14) 


where 


cohesion 


normal  stress  on  the  plane  of  failure 


p  =  excess  pressure  in  pore  water 


0  = 


angle  of  internal  friction 


(0  -  p)=  effective  stress 


Equation  14  has  two  components  [cohesion,  c,  and  friction,  (o  -  p)  tan  <?] , 
it  can  be  applied  to  all  sediments.  Shear  strength  in  sands  without  significant  amounts 
of  fine  silt  and  clay  is  defined  by  the  friction  component  (r  c  ,  these  are  "cohesion- 
less”  sediments).  Most  silt-clay  sediments  have  both  cohesion  and  friction  (under 
normal  stress).  A  few  clays  may  have  no  angle  of  internal  friction,  in  which  case  the 
shear  strength  is  defined  by  cohesion  alone. 

Equation  14  is  partly  empirical  m  that  the  exact  mechanic:  of  failure  are  not 
completely  understood,  and  it  is  not  aiw  >ys  possible  to  separate  the  contributions  of 
the  two  components  to  ultimate  shear  strength.  However,  Schmertmann  and  Osterberg 
(1960),  and  others,  have  shown  by  careful  testing  that  true  cohesion  and  friction  are 
mechanically  independent,  and  that  equation  14  is  valid  over  the  entire  range  of  strain, 
as  well  as  at  failure. 

In  sands  and  most  clays,  shear  strength  and  dy  namic  rigidity  will  increase  with 
effective  pressure,  as  indicated  in  equation  14.  The  effects  of  pressure  on  rigidity ,  how¬ 
ever,  will  not  be  considered  below,  although  very  light  effective  pressures  exist  even  in 
.,  small  laboratory  sample.  The  propagation  of  elastic  waves  in  sands  under  pressure  is 
well  understood  (Brandt,  i960,  Hardin  and  Richart,  1963;  Shell  Development  Co., 
figs.  12  and  13;  review  by  White,  1965).  Because  this  report  deals  with  surface  sedi¬ 
ments  (0  to  30  cm;  see  Part  1),  the  discussion  here  will  be  concerned  with  the  eftects  of 
cohesion  as  the  source  of  rigidity  in  high-porosity  silt-clays. 

Cohesion  is  the  resistance  to  shear  stresses  which  can  be  mobilized  between 
adjacent,  fine  particles  which  stick,  or  coheie,  to  each  other.  Cohesion  is  considered  to 
be  an  inherent  property  of  fine-grained,  clayey  sediments  which  is  independent  ot 
stress,  it  is  caused  by  physio-chemical  forces  of  an  interparticle,  intermolecuiar,  and 
intergranular  nature.  Some  important  components  and  aspects  of  cohesion  which  affect 
ngidity  are  as  follows. 


1.  Clay  particles  are  surrounded  by  layers  of  adso  bed  w.'.et  through  which  they 
interact  with  other  particles.  The  amount  of  pore  water,  the  distance  between  particles, 
and  the  number  of  interparticle  contacts  are  important  influences  on  cohesion. 

2.  Interparticle  forces  of  the  London-van  der  Waal  and  Coulombic  type,  the 
positive  and  negative  charges  on  the  faces  and  edges  of  clay  particles,  and  the  type  of 
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ions  adsorbed  on  the  clay  surfaces  and  in  the  diffuse  ion  layer  in  the  adjacent  adsorbed 
water  arc  important  contributors  to  cohesion. 

3.  At  points  of  near-contact  between  clay  particles  there  is  often  bonding  of  the 
nature  of  cementation,  especially  in  the  presence  of  iron  oxides,  calcium,  silica,  and 
other  minerals  in  solution  in  in.erstitial  waters.  Where  sediments  have  been  exposed 
to  overburden  pressures  there  is  apt  to  be  pressure-pomt  solution  and  redeposition. 

4.  The  structure  of  the  mass  of  clay  particles  is  important;  for  example,  it  has 
been  demonstrated  that  the  flocculated,  or  “cardhouse,"  structure  (Part  I,  fig.  Id)  is 
the  strongest;  these  structures  are  largely  determined  by  interpartidc  forces  and  the 
number  of  interparticle  contacts. 

5.  Differing  clay  minerals  affect  cohesion  because  of  particle  size  and  differing 
interparticle  forces;  for  example,  Na-montmorillonitc  has  stronger  cohesive  bonds  than 
kaolinitc. 

6.  Shear  stress  in  clayey  sediments  occurs  between  particles  and  not  through 
them,  near-contact  points  will  deform  elastically,  or  plastically  (depending  on  stress), 
by  an  amount  sufficient  to  sustain  the  effective  stress. 

7.  The  rate  ol  deposition,  or  age,  of  a  clayey  deposit  is  an  important  factor  in 
cohesive  strength,  li  has  been  shown  that  when  the  interparticle  bonds  of  a  slowly 
deposited  sediment  are  broken,  only  a  part  of  the  original  strength  is  regained  (thixo¬ 
tropic  regain),  especially  when  the  original  strength  is  considered  to  be  of  the  nature 
of  interparticle  cementation. 

8.  Cohesion  is  mobilized  at  very  small  strains  relative  to  the  frictional  component 
of  ultimate  shear  strength,  friction  may  be  negligible  at  the  strain  of  maximum  cohesion. 

9.  Cohesive  strength  decreases  with  increasing  temperature.  Lainbe  (  I960)  has 
diagrammatically  explained  the  components  of  shear  strength  (fig.  14). 


figure  14.  Components  of  shear  strength  of  sediments 
(from  Lam  be,  I960).  See  text  for  discussion. 

The  value  reported  as  the  shear  strength  of  a  tested  sediment  is  the  highest 
point  on  the  upper  stress-strain  curve  ("combined,  as  measured”  in  fig.  14).  This 
value  is  the  additive  result  of  the  various  components  of  shear  strength.  Cohesion  is 
mobilized  at  very  small  strains,  after  it  is  destroyed  it  ceases  to  contribute  to  overall 
shear  strength.  In  this  connection,  some  investigators  (e.g.,  Schmertinann  and  Osterberg, 
I960,  Seed  et  al.,  I960)  indicate  a  more  gradual  decline  in  the  ciuve  for  cohesion, 
which  is  logical  because  increasing  pressures  force  particles  closer  together,  thus  increasing 
some  interparticle  forces.  Dilatancy,  which  results  from  particle  interference,  causes  a 


tendency  to  volume  increase  and  more  shear  force  to  overcome.  After  additional 
strain,  there  is  no  further  tendency  toward  volume  increase;  thus,  the  dilatancy  com¬ 
ponent  is  overcome,  and  the  important  component  is  friction.  When  the  overall 
stress-strain  curve  becomes  horizontal,  the  only  component  of  shear  strength  is  friction, 
including  particle  interference,  which  lessens  as  the  clay  platelets  tend  to  align  them¬ 
selves  with  their  long  axis  parallel  to  the  direction  of  shear. 

To  destroy  the  sediment  structure,  shear  stress  must  first  break  down  the  com¬ 
plex  of  interparticle  forces  and  cementing  bonds  outlined  above.  When  a  clay  sediment 
structure  and  cohesion  are  destroyed  by  stress,  porosity  is  reduced  under  additional 
pressure,  so  that  further  tests  of  sound  velocity  are  not  realistic  and  results  should  not 
be  applied  to  real  sediments.  This  is  why  velocity  measurements  of  silt-clay  sediment 
samples  under  laboratory  pressures  (as  in  the  consolidation  test)  cannot  be  directly  related 
to  the  same  pressure  levels  in  a  natural  sediment.  Conversely,  any  artificial  clay-silt 
sediment  composed  in  the  laboratory  will  not  have  the  interparticle  bonds  (cohesion)  of  a 
natural  sediment,  especially  those  sediments  of  the  deep  sea  which  have  accumulated 
slowly  over  geologic  time. 

A  recent  study  with  artificially  sedunented  clays  dramatically  demonstrated 
some  of  the  effects  of  sediment  structure,  interparticle  bonding,  thixotropic  regain  of 
strength,  and  mineralogy  on  dynamic  rigidity  (Cohen,  1 968).  Some  other  important 
aspects  of  this  study  were  discussed  in  the  section  “Rocks  and  Sediments  as  Elastic 
Media."  Cohen  used  a  torsional  oscillator  resonance  technique  to  measure  complex 
dynamic  rigidities  in  mixtures  of  kaolinite  and  bentonite  in  distilled  water  at  various 
frequencies.  His  results  included  the  following 

1.  Kaolinite 

a.  Kaolinite  concentrations  of  about  32  percent  (by  weight!  in  distilled  water 
formed  flocculated  stiucturcs  with  densities  and  porosities  comparable  to 
high  porosity  sediments.  The  elastic  portion  of  complex  rigidity ,  /a,  in  one 
experiment  increased  from  0.8  dynes, cm-  X  105  alter  one  day,  to 

2.16  dynes/cm-  X  105  after  5  days. 

b.  When  Calgon,  a  deflocculating  agent,  was  added  to  the  above  mixture,  the 
structure  became  dispersed,  no  elastic  rigidity  was  measured,  and  the  mix¬ 
ture  behaved  as  a  Newtonian  fluid. 

c.  When  35.5  ppt  of  NaCI  was  added  to  the  Calgon  mixture,  the  clay  floccu¬ 
lated  and  rigidities  were  about  the  same  as  before. 

2.  Bentonite 

a.  Bentonite  concentrations  of  about  19  percent  in  distilled  water  formed 
flocculated  structures  in  which  densities  were  lower  and  porosities 
higher  than  in  natural  sediments.  Elastic  rigidities  increased  in  a  typical 
sample  from  1.33  dynes/cm2  X  105  after  one  day.  to  4.ts7  dynes/cm2  X 
I05  after  3  days  (with  little  or  no  increase  in  the  next  few  days). 

Cohen’s  results  demonstrated: 

1 .  the  quantitative  rigidities  mob.lized  in  flocculated  clay s,  after  a  few  days,  as 
the  result  of  interparticle  bonding  and  thixotropic  regain  of  strength  (after  defloccula- 
tion).  and  the  increase  of  rigidity  with  time, 

2.  the  destruction  of  interparticle  bonding  by  addition  of  a  deflocculating  agent, 
and  resumption  of  bonding  in  the  presence  ol  an  electrolyte;  and  that  flocculated 


structures  of  this  type  have  measurable  dynamic  rigidities,  and  dispersed  structures 
behave  as  fluids  or  suspensions, 

3.  that  bentonite  (dominantly  montmorillonite)  forms  stronger  structures  than 
docs  kaolinite,  a  well-known  f3ct  in  soil-mechanics  research  ( c.g Warkentin  and  Yong, 
1962;  review  by  Me’ade,  1964). 

Slow  rates  of  deposition  and  great  age  of  deposits  cause  increases  in  sediment 
structural  strength  due  to  increased  intergrain  bonding  having  the  nature  of  cementation 
(Leonards  and  Ramiah,  l959;Bjerrum  and  Wu,  1960,  Bjerrum  and  Lo,  1963;  Meade, 

1963  ,  Leonards  and  Alischaeffl,  1964)  Such  increase  in  strength  is  reflected  in  the  lack 
of  appreciable  reduction  in  porosity  with  overburden  pressure  in  deep-water  silt-clays 
(rev'ewed  in  Part  I),  and  in  the  relatively  high  shear  strengths  in  these  sediments 
(Moore,  1961,  1964;  Hamilton.  1964,  Richards  and  Hamilton,  1967).  The  maximum 
rigidities  measured  by  Cohen  (1968)  were  about  2  d>  nes/cm2  X  106.  In  the  San  Diego 
Trough,  m  situ  determinations  of  shear-wave  velocities  and  densities  resulted  in  com¬ 
puted  least  values  of  about  1  dyne/ern-  X  10R  (Hamilton  ft  a!.,  1969).  In  deep-sea 
clav«,  computed  values  of  rigidity  were  about  5  Jyncs/em2  X  10®.  The  increased  values 
in  the  natural  sediments  probably  include  additional  rigidity,  resulting  from  the  effects 
of  age  and  bonding,  which  has  the  nature  of  cementation.  Schreiber  ( 1968)  has  demon¬ 
strated  the  effects  of  cementation  in  volcanx  glass-ash  lay  ers  in  a  natural  marine  sediment 
from  the  deep  Caribbean  Sea.  These  layers  had  unusually  high  compressionaj  velocities 
and  sediment  strength  (cohesion).  The  addition  of  HC1,  or  loss  of  water  by  drying, 
resulted  in  marked  loss  of  strength. 

Hardin  and  Richart  (1963)  have  demonstrated  that  in  sands,  shear-wave  velocities 
are  independent  of  grain  size  and  size  distribution,  except  as  they  affect  porosity ;  and 
that  dynamic  rigidity  decreases  with  increasing  porosity.  However,  sands  composed  of 
angular  grains  have  higher  rigidities  (at  the  same  porosities)  than  do  those  composed  of 
round  grains,  apparently  because  of  increased  resistance  of  the  angular  grains  to  inter¬ 
grain  movements. 

Shumway  (1960)  related  the  amount  of  sand  (fraction  of  a  unit  volume)  in  his 
samples  to  rigidity  under  the  assumption  that  these  grains  were  in  contact  and  that  the 
greater  their  volume,  the  greater  the  elastic  rigidity  of  the  sediment  structure.  This  is  an 
invalid  assumption;  in  the  higher  porosity  sediments  (even  if  ptesent  in  significant 
amounts)  such  larger  grains  are  not  in  contact  (Part  I,  fig.  1 0  as  pointed  out  (in  the  con¬ 
text  of  elasticity)  by  Kozlov  (1962).  Shumway  also  assumed  that  rigidity  was  zero  at 
80  percent  porosity,  which  is  incorrect  in  the  laboratory  (Cohen,  1968)  or  in  situ 
Hamilton  et  ai,  1969;  discussions  above). 

QUANTITATIVE  EFFECTS  OF  RIGiDITY  ON  COMPRESSICNAL  W  AVE 

VELOCITY 

The  effects  of  rigidity  on  values  of  comprcssional-wave  velocity  in  natural 
marine  sediments  are  small,  but  significant.  Computations  were  made  which  involved 
equation  I ,  measured  values  of  cc-mpressional-wave  velocities  and  densities,  and  either 


values  of  p  computed  from  equation  6,  when  shear-wave  velocities  were  known,  or 
those  computed  aftei  an  initial  computation  of  the  bulk  modulus.  The  results  were: 

Sediment  Type  f'p  (m/sec)  Remarks 


with  p 

without  p 

Fine  sand  (laboratory) 

174(1* 

1 7  Ih 

Shell  Development  Co.  Jaia 

Medium  sand  (in  sun ) 

1798’ 

1784 

Hamilton  et  al.  ( 1969) 

Fine  sand  (corn  terrace) 

1742” 

1696 

Silty  clay  (abyssal  hills) 

1507” 

1486 

*  L  1  s,  p  measured,  p 

nnipulcd 

*T  ,p  measured,  k  ,  and  then  p  computed 


In  general,  the  'omputations  indicate  that  the  presence  of  dynamic  rigidity  in 
sands  is  apt  to  raise  comptesstonal-wave  velocities  on  the  order  of  I  to  3  percent;  in  deep- 
sea  clays  the  average  increases  should  be  about  1  to  2  percent 


CONCLUSIONS  REGARDING  RIGIDITY 

The  data  presented  in  previous  sections  indicate  that  almost  all  marine  sedi¬ 
ments  which  have  mineral-to-mineral  contacts  (as  in  sands),  or  flocculated  clay  structures, 
have  small  but  significant  rigidities  which  allow  transmission  of  shear  waves.  This  is  true 
even  in  artificial  clays  in  distilled  water  (Cohen,  1968).  Thus  an)  equation  for  elastic- 
wave  propagation  in  natural  marine  sediments  which  does  ot  provide  for  rigidity  (as  the 
W  ood  equation,  5)  should  be  abandoned. 

It  is  apparent  from  the  discussion  of  cohesion  that  dynamic  rigidity  cannot  be 
computed  for  a  system  such  as  deep-sea  clay  from  us  physical  components  (r  e.,  given 
porosity,  mineralogy,  and  pore-fluid  composition)  One  might  be  able  to  compute  some 
ol  the  interparticle  forces  using  clay -mineral  technology,  but  the  bonding  resulting  fiont 
age  of  the  deposit  and  “cementing”  effects  could  not  be  computed.  However,  enough 
information  is  at  hand  to  reasonably  predict  values  of  dynamic  rigidity ,  given  the  sedi¬ 
ment  type  and  environment,  plus  physical  prupci  lies  such  as  density,  porosity,  and 
mineralogy. 

Although  there  are  "usable”  empirical  relationships  between  rigidity  and  other 
physical  properties  ( figs  1  5,  lb),  a  belter  procedure  for  deriving  a  value  for  rigidity  (lack¬ 
ing  shear-wave  velocity)  is  to  use  equation  8  and  values  for  density,  compressions!  veloc¬ 
ity,  and  a  computed  value  fot  the  bulk  modulus  especially  for  the  continental-terrace 
environment  where  scatter  is  great  between  rigidity  and  any  common  physical  properly. 
Lacking  values  for  any  physical  properties,  the  average  values  of  rigidity  for  the  esti 
mated  sediment  type  in  the  particular  environment  should  be  used.  If  grain-size  data, 
density,  or  porosity  are  known,  missing  values  of  density,  porosity,  and  velocity  can  be 
determined  for  23°C  and  1  atmosphere  by  using  procedures  discussed  in  Far!  I,  in  situ 
values  should  be  computed  according  to  the  procedures  of  Fart  III 

Empirically,  the  best  index  to  rigidity  is  p  I  in  both  the  continental  terrace 
and  abyssal-plain  environments  (fig.  16).  In  the  abyssal-hill  environment,  the  best  index 
is  percent  clay  size  (or  sand  plus  Sill;  see  fig.  I  5).  Regression  equations  aie  included  in 
Appendix  C  for  these  three  relationships 


POISSON'S  RATIO 


When  rigidity  is  zero,  no  shear  wave  can  be  transmitted,  and  Poisson's  ratio  is 
that  of  a  fluid  ot  suspension,  0.50  (eq  10)  As  discussed  in  the  preceding  section,  most 
natural  marine  sediments  possess  rigidity  and  transmit  shear  waves;  therefore,  most  sedi¬ 
ments  have  values  of  Poisson’s  ratio  less  than  0.50. 

As  noted  previously,  many  literature  values  of  compressional-wave  velocities  in 
"saturated"  sediments  (especially  sands)  must  be  disregarded  because  the  measurements 
were,  apparently,  made  in  sediments  at  less  than  lull  saturation,  where  shear-wave  veloci¬ 
ties  are  slightly  higher,  bur  comprcssional  velocities  arc  far  too  low.  As  a  result,  many 
published  values  for  Poisson’s  ratio,  especially  in  soil  mechanics  literature,  are  too  low 
because  the  ratio,  f'  /fj.  determines  Poisson's  ratio  (eq.  10). 

In  the  references  tu  table  3  there  are  1  1  cases  where  comprcssional-  and  shear- 
wave  velocities  in  the  laboratory  or  in  situ  allow  computations  of  Poisson’s  ratio;  these 
are  within  the  range  of  this  ratio  as  computed  in  the  present  study  (tables  1  and  2). 

Average  values  and  standard  errors  of  Poisson’s  ratio  are  listed  in  tables  I  and  2. 
Maximum  and  minimum  values  (not  listed)  indicate  the  following  ranges  and  averages  in 
some  principal  sediment  types  and  environments. 


Environment 

Sediment 

Poisson’s  Ratio 

No.  of 

Samples 

Max. 

Min. 

Avg. 

Continental  Terrace: 

Sands  (all  grades) 

0496 

0416 

0.470 

13 

Clayey  sdt 

0.499 

0.447 

0.478 

40 

Abyssal  Plain  (Tutbidilc) 

Sdt-clays 

0  496 

0.466 

0.484 

52 

Abyssal  Hill  (Pelagic): 

Silt-clays 

0.499 

0.467 

0.487 
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Tire  data  in  the  table  above  are  for  laboratory  conditions.  Anderson  and 
Sehreiber  ( 1 965 ).  Anderson  and  Liebermann  (I96X),  and  Nate  and  Drake  ( 196?)  have 
noted  that  Poisson's  ratio  does  not  change  substantially  with  pressure  and  temperature  in 
rocks.  Using  the  methods  discussed  in  Part  III.  computations  comparing  laboratory  and 
in  situ  values  of  Poisson’s  ratio  for  a  given  high-porosity  deep-sea  silty  clay  indicate  this  is 
also  true  of  marine  sediments.  W  ith  temperature  and  pressure  the  only  variables, 
Poisson’s  ratio  in  the  laboratory  at  23°C  and  1  atmosphere  is  0.486;  at  a  water  depth  of 
6000  m  (pressure:  626  kg  cm-  ;  temperature:  1 .5°C),  Poisson’s  ratio  is  0.490. 

Velocity-Frequency  Relationship 

The  subject  of  velocity  dispersion  (dependence  of  velocity  on  frequency  )  is 
important  in  geophysics  and  underwater  acoustics.  Questions  concerning  basic  theories 
and  models  for  propagation  of  elastic  waves  in  various  media  cannot  be  resolved  until 
various  questions  are  answered:  does  velocity  dispersion  exist;  and,  if  so.  to  what  extent, 
in  which  media,  and  over  what  frequency  range'1  The  subject  of  velocity  dispersion  in 
satuiated  sediments  was  recently  raised,  again,  by  Hampton  (1967)  who  reported  disper¬ 
sion  in  laboratory  measurements  in  aruficia!  sediments. 
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Because  well-known  and  valid  theoretical  models  exist  on  both  sides  of  the 
question  of  dispersion  m  earth  materials,  the  answer  must  lie  in  experimental  work  which 
will  set  parameters  for  these  models.  It  is  the  opinion  of  the  author  that  the  expert 
mental  evidence  strongly  indicates  that  in  the  frequency  range  from  a  lew  hertz  to  the 
megahertz  range,  there  is  no  measurable  dependence  of  velocity  on  frequency  in  earth 
materials,  including  almost  all  natural  sediments  and  rocks. 

In  any  discussion  of  the  evidence  concerning  velocity  dispersion,  a  cleat  distinc¬ 
tion  should  be  made  between  various  media.  The  broad  categories  are: 

1.  rocks  with  little  or  no  porosity. 

2.  porous,  saturated,  and  unsaturaied  rocks, 

3.  fully -saturated  sands  with  mmeral-to-imneral  gram  contacts, 

4.  higher-porosity,  fully  saturated  sediments  with  a  structure  possessing  cohesion, 

5.  partially  saturated  sediments  (natural  or  laboratory )  containing  gas  or  air  in 
pore  spaces,  and 

6.  clay -silt  “sediments"  which  arc  suspensions,  especially  those  formed  artificially 
in  the  laboratory  with  a  deflocculating  agent. 

The  experimental  evidence  or  velocity  dispersion  falls  into  the  categories  of 
work  in  the  laboratory  on  rocks  and  on  natural  and  artificial  sediments,  and  in  situ  on 
rocks  and  sediments  Measurements  of  both  velocity  and  sound  absorption  are  involved 
The  references  cited  below  are  recent  examples  of  this  experimental  work,  no  compre¬ 
hensive  review  is  intended. 

A  number  of  investigators  have  measured  compressional-  and  shear-wave  veloc¬ 
ity  and  absorption  in  rocks  (laboratory  and  in  situ).  Their  common  conclusion  is  that 
there  is  no  (or  negligible)  measurable  velocity  dispersion  in  the  range  from  seismic  fre¬ 
quencies  into  the  megahertz  range,  and  that  the  specific  dissipation  constant,  1  IQ,  is  in¬ 
dependent  of  frequency.  Examples  include  work  by  Wyllie  era/.  (195b),  Birch  (1961), 
Peselnick  and  Outerbridge  (1961),  McDonaJ  et  ol.  (1958),  White  (1965),  Press  (1966). 
Bradley  and  Fort  (1966)  have  good  resumes  of  much  of  the  evidence. 

No  velocity  dispersion  was  measured  in  artificial  sands  in  the  laboratory  by 
Hardin  and  Richart  (1963),  Nolle  et  al.  (1963).  and  Schorl  (1963);  these  measurements 
included  a  frequency  range  from  200  Hz  to  1  MHz. 

In  soil-mechanics  investigations  in  siru  in  sands,  low-frequency  vibrations  were 
used  in  studies  by  Barkan  ( 1962)  and  Jones  ( 1958)  to  measure  shear-wave  velocities 
(Jones'  measurements  also  included  clay-silt ).  No  velocity  dispersion  was  measured  in 
the  frequency  range  from  10  to  400  Hz. 

Compressional-velocity  studies  in  high-porosity,  deep-sea  clay  in  the  Nor'h 
Pacific  by  the  writer  at  200  kHz  (laboratory),  by  Schreiber  (1968b)  and  Horn  et  al. 
(1968b)  at  400  kHz  (laboratory),  and  by  Fry  and  Raitt  (1961;  seismic  measurements  at 
sea)  are  all  in  reasonable  accord  (less  than  1  percent  differences).  The  average  of  com- 
pressional  velocity  for  abyssal  hill,  silty  clay  (Part  I.  table  2)  Is  I  507  m/sec.  An  average 
velocity  of  the  top  of  10  cores  off  Hawaii  (Schreiber,  1968a),  in  the  same  material,  is 
1504  m/sec.  The  average  ratio,  velocity  in  sediment/velocity  in  water,  at  200  kHz  (Part 
I,  table  2)  for  abyssal-hill  silty  clay  is  0.985  (max  1 .006,  min  0.973).  The  ratio  comput¬ 
ed  by  Fry  and  Raitt  ( 1961 ,  table  1 )  for  this  sediment  type  was  0.974  (see  discussion 
in  Part  I). 

In  a  continuing  program  to  measure  compressional  velocity  and  attenuation  in 
siw.  the  writer  measured  no  velocity  dispersion  in  shallow-water  sands  (three  stations), 
or  sill-clay  tone  station)  at  14,  7,  and  3.5  kHz. 

In  a  laboratory  study  of  complex  rigidity  {n  +  /)/)  in  artificially  sedimented 
clays,  Cohen  ( 1 968)  demonstrated  that  both  components  of  complex  rigidity,  ju  and  rp  , 
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were  independent  of  frequency  in  the  range  8.6  to  43.2  kHz.  When  a  deflocculating 
agent  W3s  added,  the  flocculated  structure  was  dispersed,  the  material  lost  ail  rigidity,  p, 
and  the  mixture  behaved  as  a  Newtonian  viscous  fluid  in  winch  i; 1  was  linearly  depend¬ 
ent  on  frequency.  The  addition  of  35.5  ppt  of  NaCI  caused  reflocculation  and  both  ele¬ 
ments  -if  complex  rigidity  were  the  same  as  before. 

There  are  several  important  conclusions  to  be  drawn  from  Cohen's  study  (in 
addition  to  those  discussed  in  previous  sections): 

I  Suspensions  (without  structure)  do  not  behave  as  flocculated  day  structures, 
and  almost  all  natural,  high -porosity  silt-clays  have  lilts  type  of  structure  fsec  section  on 
“Causes  of  Rigidity”). 

2.  If  the  clastic  portion  of  complex  rigidity. /a.  is  independent  of  frequency,  there 
is  no  shear  velocity  dispersion  in  the  flocculated  clay. 

Hampton  ( 1 967),  in  a  laboratory  study  of  artificial  sedunents,  reported  a  4  to  6 
percent  increase  in  compressional  velocity  from  3  to  2(X>  kHz  (this  amounts  to  60  to  90 
m/scc  in  silt-clays);  his  ratios,  velocity  in  sediment/velocity  in  water  ( 1967,  fig  I  1),  indi¬ 
cated  sediment  velocities  as  low  as  7  percent  below  that  in  water  at  lower  frequencies 
Hampton  referred  to  Ament  ( 1 953 )  for  theoretical  support  of  his  conclusions. 

Ament’s  (1953)  theoretical  approach  to  sound  propagation  in  gross  mixtures 
involved  a  true,  viscous  suspension  in  which  particles  were  not  m  contact  and  in  which 
permeability,  viscosity,  and  scattering  of  sound  by  the  particles  were  involved.  He  com¬ 
pared  his  approach  with  experiments  by  Urick,  who  measured  compressional  velocities 
through  suspensions  of  mercury,  biomofonn,  and  kaolinite  (with  a  deflocculating  agent) 
in  water.  Sutton  et  at  (1957)  and  Laughton  ( 1957)  discussed  Ament’s  equations  and 
concluded  that  they  could  not  be  applied  to  natural  marine  sediments,  the  author  agrees 
with  this  conclusion.  As  discussed  above,  subsequent  measurements  of  elastic-wave 
velocities  in  the  laboratory  and  field  do  not  support  any  theory  involving  significant 
velocity  dispersion  from  low  (seismic)  frequencies  to  the  megahertz  range. 

Hampton’s  values  of  the  ratio,  velocity  in  sediment  to  velocity  in  water  ( 1967,  fig, 

1 1 )  not  only  varied  with  frequency  but  were  anonymously  low  fat  lower  frequencies). 

His  values  for  artificial  kaolinite  "sediments,"  for  example,  in  the  fiequency  range  8  to 
40  kHz  (0.93  to  0.97),  can  be  Compared  (at  similar  concentrations)  to  those  of  Cohen 
fl968)  in  kaolinite  suspensions  and  flocculated  clay  in  the  same  frequency  range  (0.97  to 
0.99),  and  those  of  Urtck  ( 1947)  in  a  kaolinite  suspension  at  I  Mllz  (0.97  to  0.99). 
Shumwav  ( 1 958,  1 960)  measured  a  ratio  of  0.97  in  a  red-clay  slurry.  An  average  utio 
for  six  of  Sit  uni  way’s  (I960)  high  -porosity  sediments  in  the  San  Diego  Trough  (20  to  40 
kllz I  waa  0.98.  Im  equation  7  ( i 960,  fie.  i ,  p  io i)  predicts  an  average  least  value  ot 
0.97  for  his  sediments.  The  author  ( 1956)  measured,  in  situ,  in  the  Sail  Diego  Trough, 
average  values  of  0.98  at  100  kH2  Fry  and  Raitl  ( 1967)  computed  an  average  valu.  of 
0.974  for  deep  Pacific  surface  sediments  at  seismic  frequencies.  In  the  Atlantic,  Houtz 
and  Ewing  (1964 )  used  the  same  measuring  lechmques  thai  were  used  by  Fry  and  Raitt 
to  obtain  sediment  surface  velocities.  At  three  stations  (Stas.  4,  5,  and  8)  favored  in 
later  discussions,  Houtz  and  Ewing  measured  values  of  1517.  1 532,  and  1517  m/sec, 
which  when  divided  by  appropriate  values  for  bottom-water  velocity  yield  an  average 
ratio  of  0.98.  Table  2  in  Part  1  lists  average  ratio  values  in  deep-water  silt-clays  (at  200 
kHz)  in  the  range  0.98  to  0.99.  In  summary,  the  evidence  from  laboratory  and  field, 
over  a  wide  range  of  frequencies,  does  not  support  Hampton’s  low  ratio  values 

When  air  or  gas  bubbles  (from  decaying  organic  material)  are  trapped  in  pore 
spaces  within  any  sediment  they  have  a  marked  effect  on  both  compressional  velocity 
and  attenuation,  depending  on  the  concentration  and  size  ol  bubbles,  and  the  frequency 
(r.£  ,  Meyer,  1957):  velocity  vanes  with  frequency  (velocity  is  usually  too  low)  and 
attenuation  is  apt  to  be  high.  Ad  though  Hampton  recognized  these  facts,  and  attempted 


u>  remove  gas  an J  air  bubbles  I  tom  Ins  aitilicial  sediments.  Ins  data  indicate  the  proba¬ 
bility  ot  an  jnd.-oi  gas  bubbles  within  lus  materials. 


Summary  of  Factors  Affecting  Compressional- 
Wave  Velocity 

1  o  summarize  the  ellects  ol  many  complex .  variable  (actors  on  compressional 
velocity  .  it  is  instructive  to  separate  equation  I  into  two  components  and  to  expand 
density  .  p,  as  pei  equation  4.  and  k  as  pci  equation  Ida. 
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The  variable  factors  contributing  to  compressional  wave  velocity  are  thus 
Bulk  moduli  of  pote  water  ).  mineials  (xs).  and  frame  ). 

Porosity  (ul.  in  a  unit  volume,  porosity  equals  the  volume  of  pote  water 
and  ( I  -ii )  equals  the  volume,  or  concentration,  of  mineral  grams. 

Densuy  of  pore  water  (pw. ),  and  minerals  (psl 
Rigidity  modulus  (p) 

Given  laboratory  conditions  of  constant  tempetalure  and  pressure,  and  assuming 
that  pore  water  is  the  same  (true  within  narrow  limits;  see  Appendix  B).  the  effects  of 
changes  in  porosity  and  mineralogy  are  t'ms  (summarizing  previous  discussions): 

1.  Changes  in  mineralogy  only  cause  changes  in 

a.  the  bulk  modulus,  k,  through  changes  in 
( 1  )  the  bulk  modulus  of  minerals.  ks 
(2)  the  frame-bulk  modulus,  * 

b.  bulk  density, p,  through.  (1  -u)ps 

c.  the  rigidity  modulus,  n\  in  sediments  hat  •  cohesion,  some  minerals  have 
siiongei  interpariieio  bonds  (see  "Causes  ol  Rigidity") 

2.  Changes  in  porosity  only  cause  changes  in 

a.  the  bulk  modulus,  k  ,  through  changes  in 
( I  )  the  I  rame-bulk  modulus,  k  ,  ( fig  2 ) 

J 

(2)  the  denominator  of (d(eq.  13a):  n  (*s  -k  j 

b.  bulk  density ,  p.  through  the  products  (eq.  4):  npw  and(l-rrlp 

c.  the  rigidity  modulus,  ju.  because  higher  porosities  lead  to  fewer  tnterpaiticle 
contacts  in  sands  and  in  si) l -clay  s  (having  cohesion,  see  "Causes  of  Rigidity") 

Die  causes  of  variations  in  rigidity  were  discussed  in  a  previous  section,  but  the 
results  of  these  variations  can  be  readily  seen  when  equation  I  is  separated  into  its  two 
components,  as  above.  When  rigidity,  p,  is  zero  or  negligibly  small,  the  component 

4/3  P 
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ilni|)v  mil  and  lln'  icniamdci  of  equation  I  is  (In'  equation  I'm  eomprossioi,:il  velocity  ill  a 
liquid,  or  any  medium  without  rigidity  It  is  the  same  as  the  Wood  equation  ( 5 )  it  the 
medium  is  a  true  suspension  with  no  sediment  sltuetute.  in  wlneli  ease  (cq.  3b)  can 
be  used  for  the  bulk  modulus,  k.  If  there  is  a  sediment  struetute,  x  should  include  a  frame- 
bulk  modulus,  hj.  Because  almost  all  marine  sediments  have  a  definite  structure  and 
possess  rigidity .  the  Wood  equation,  which  served  to  approximately  dc: line  cotnprcssional 
velocity  in  much  catly  work,  should  be  abandoned  in  favor  of  equation  1 ,  the  classic 
equation  for  compressions!  velocity. 

1  he  interaction  of  all  of  the  abuse  actors  r-  sults  in  almost  all  higher-porosity 
sediments  having  eoniprcssional  velocities  less  than  in  water  (poic  water  or  bottom  water 
above  the  sea  floor)  unless  rigidity  is  unusually  high  This  low  velocity  phenomenon  is 
mostly  due  to  the  low  rigidities  in  high-poiosily  sediments,  and  the  dominant  ellccts  ot 
high  wjtct  compressibilities  (oi  loss  incompressibilities)  relative  to  mineral  moduli  (see 
Appendix  If  for  typical  values) 

Summary  Statement  Regarding  Elasticity  of 
Marine  Sediments 

To  prove  that  a  medium  responds  as  an  clastic  body,  ii  is  necessary  to  use  the 
theory  of  elasticity  and  its  equations  (eqs.  1 ,  2,  5-1 1 }  to  compute  unmeasured  elastic  con¬ 
stants  (without  empirical  constants  or  assumptions)  and  verily  these  by  measurements. 
Further,  values  for  elastic  constants  should  be  computable  from  the  basic  components  of 
the  medium,  and  static  and  dynamic  values  should  be  theoretically  related.  All  of  these 
requirements  have  been  satisfied  in  the  ease  of  rocks  ( c.g. ,  Birch,  1961  .  Simmons  and 
Brace,  1 965 ,  Anderson  and  Liebcnnann,  1968).  It  is  a  conclusion  of  this  report  that 
these  requirements  can  also  be  satisfied  in  the  case  of  saturated,  porous  sediments 

Density  and  compressional-wavc  velocity  were  measured  m  the  present  study  A 
computed  value  for  the  bulk  modulus  was  used  as  the  third  constant  required  to  compute 
the  other  elastic  constants  with  the  equations  of  elasticity.  This  bulk  modulus  was  com¬ 
puted  from  the  porosity  and  the  bulk  moduli  of  the  components  of  the  medium  (pore 
tlu id.  minerals,  and  a  frame-bulk  modulus)  following  (j.issmaiui’s  equations  ( 1 3a,  Appen¬ 
dix  A).  This  method  is  thus  based  on  theory  without  empirical  factors  or  constants. 

The  computations  of  elastic  constants  as  discussed  in  this  report,  compared  with 
other  laboratory  and  in  situ  measurements  (such  as  shear  velocity,  tabic  3),  indicate  that 
the  equations  of  elasticity  can  be  used  to  derive  reasonable  values  for  those  elastic  con¬ 
stants  not  measured  If  so,  the  computed  values  in  tables  1  and  2  predict  these  constants 
for  ihe  major  marine-sediment  types. 


REVERSE  SIDE  BLANK 
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APPENDIX  A: 

DERIVATIONS  OF  GASSMANN’S  EQUATIONS 

The  derivation  of  Gassmann’s  equation  t  J ‘>5  1  ,  p  15,  paragraph  No.  59) 

K'/  +  y 

K  -  K  -  Q  - - '■ — 

1  «.  *  Q  n  (ks  -  kk.) 

is  indicated,  hut  not  completely  shown  in  Ins  original  (German)  publication  The  follow¬ 
ing  derivation  (by  the  author)  is  completed  because  ot  its  importance  to  the  computa¬ 
tions  of  this  repot t.  and  to  make  it  more  readily  available  to  interested  readers 

Gassmann’s  equation,  above,  is  not  empirical;  it  follows  from  the  basic  assump¬ 
tions  that  a  unit  volume  of  a  saturated  porous  medium  responds  elastically  under  the 
light  pressures  of  a  sound  wave  In  this  “closed  system"  no  pore  water  leaves  or  is  added 
to  the  unit  volume  during  passage  ol  the  sound  wave. 

The  notations  of  Cassmann  are  changed  as  follows 


Item 

Gassmann 

Changed  To 

Total  pressure 

P 

No  change 

Hydiostatic  pressure 

P 

Effective  pressure,  or  pressure  on  frame 

P 

P 

e 

Bulk  modulus  of  system 

K 

No  change 

Bulk  modulus  of  pore  watei 

K 

V 

Bulk  modulus  of  mineral  solids 

A 

K 

Kc 

Bulk  modulus  of  mineral  frame 

h 

7 

Total  volume 

T 

No  change 

V  olume  of  pore  water 

T 

T 

Volume  of  minca!  solids 

T 

r 

s 

I .  With  these  changes  in  notation,  peitmenl  equations  and  "running"  translations 
from  Gassmann  ( 195  1 .  p  I  5)  are 

"(5b)  In  order  to  determine  r. .  the  unit  volume.  I  .  will  be  exposed  to 
an  additional  total  (hydrosiaticj  pressure.  A P.  which  is  effective  on  all 
surfaces  .  A/5  is  icsclvcd  ml-;  two  cnn'ipui'icnr-s  A/y  is  that  paii  oi 
AP  which  is  effective  on  ah  surfaces  of  the  volume  I  ,  and  determines 
the  addi'mnul  hydrostatic  pressure  in  the  ’ores  A P ,  is  the  remainder 
of  A P  which  is  effective  on  those  portions  of  the  surfaces  occupied  by 
the  mineral  solids  of  the  frame 

a  A/’  =  IP,  +  A  P  (A- la) 

I  he  desired  bulk  modulus,  k  .  of  the  closed  system  is  defined  bv  the 
equation 

b 


AT 

V 


A/' 


AT 


a  r 


(A-ib) 


( A  - 1  c ) 
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describes  the  compression  ot'  the  pore  walet  (or.  Itoni  I ): 


A P,  n. k 
h  ■  »■ 


describes  the  compression  ol  the  tianie.  proceeding  troin  the  two 
components  A/1  in  (.1). 

AC  ( 1  -  ri|  A/5,  A P, 


describes  the  compression  ol  ti  e  mineral  solids  .  .  with  the  application 
of  l’,  =  (  I  /i >  I'. 


det  ines  porosity . 


I'  =  n  I 


AT  =  AT  +  Al 

s  VI 


which  indicates  that  the  system  is  closed,  therefore,  no  pore  water  can 
be  added  nor  can  any  escape.” 

Using  the  above  basic  relationships  furnished  by  (Jassmann.  the  final  equation  lot 
computing  the  system  bulk  modulus,  k .  can  be  derived  as  follows: 

2.  Equating  volume  changes 

u-  Al  A I 

A  t  __  _ v  ^ _ w 

i  ~  r  r 

from  b.  c,  and  e,  above,  and  multiplying  by  <  -  I ) 

sn  A/>.  A  P  A  P.n 

A/  .  n  e  n 

-  =  (1/0  -  +  -  +  - 


from  a  A/5  =  A P-  A P, 


reduces  to 


A/)  Il-HIA^  A  Pf,n  A/'  A  P.. 


A Ph  |'H\  -  sM.)l 


Divide  bv  A/’ 


I  _  Sv»]  +  _1_ 


K  A P  K  K 

S  H' 
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:  i 
r  1 
I  { 


)  1 


Al’„ 


K  K  AP  K  K 

S  5  W 


reduces  u> 


AI),  l"<S~ 

AP  K  K 


AP  _  K  1",KJ  “  \4>1 


AP, 


K  (K  “A  ) 
ll'  J  ' 


tquaung  pressure  .'olume-n'oduli  relationships 

a.  Compression  of  frame  (deoieases  «f  .  ..  Lq  A-l(d) 

b.  Compression  of  solids  (increases  >r) .  Lq  A-l(e) 

c.  Compression  of  water  . .  f;q.  A-l(c) 

In  the  closed  system  no  pore  water  is  added  or  escapes,  therefore 

ia  -  a b  -  ,V.  or 


*/ 


(I  ->»±Ph  ^  AP. 


yP).n 


AP,  AP  A P,  A P.  n  A P  AP,  n 

/i  +  f  //  ^  /;  _  /i 


AP  AP.  AP.  n  AP,  n 


AP,  (K  -k  )  AP,  r;  (k  ) 

<  -  J  _  ti  S  Vv 


AP  = 


AP,.n(K.  -k  IK,. 


K  (K  -K,.| 

it  s  ! 


from  eq  la.AP  =  AP  -  AP, 

1  e  /; 


AP  -  AP, 


=  -yV1(VK"-!  V 
K«-fvV 


AP  = 


AP,  // (k  -k  )k  +  AP.  k  (k  -k,)1 
/:  i  u  /  h  1  w  s  t  1 


k  Ik  -k  I 
w  s  ; 


Divide  bv  A P. 


M»  n <K  K  )  K  .  *  K  IK  “K , ) 
__  ~  V  H‘  ’  J _ w  A  / 

AP,  K  (k  -k  ) 

n  w  s  r 


(A  2) 


( A--') 
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4  liquating  equations  A-2  and  .4-3 


K 


U  -K) 


n  Ik  ~k  )  k t  +  k  ( k  ~k  J 
v  s  w'  J  V v  s  f 


let 


*V(  W  '  A 
n(lis-k„)  = 


then 


kB 


B  k  -  +  4 


a  In  -K )  /I 
w  r 

wluth,  after  dividing  bv  S,  reduces  to 


(d  to  +  V  —  \  =  K  ,.((  K  +  K  K 

\  /  *  '  **■ R!  1  1  “  s  w' 


Substitute  43  and  solve  for  k 

D 


K  K  (li.^Q) 

_ tv  S  J  ^ 

.4  +  k  -  n  +  x  (d 

J  H'  kV 


K»-Ks  V^> 


K  K  -X  K  ,  +  X  K  ,  +  X  4> 
$  w  /  w  J  w 


which  reduces  to 


which  is  Cussmann  s  c< 


K - 

f  IOC  l  «  1  c  n.».- 

sfud  tlOII  I  I  e  y  j  ,  p  .  J  ,  pul 
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APPENDIX  B: 

BULK  MODULI  OF  WATER  AND  MINERALS 

Some  values  ol  bulk  moduli  (incompressibilities)  of  seawater,  k^,  and  minerals, 
k  used  in  the  computations  in  the  main  test,  are  listed  below.  In  addition,  values  f  ir 
the  abrogate  bulk  moduli  loi  some  of  the  more  important  types  of  sediment  mineral 
aggregates  are  listed.  These  were  computed  by  the  Voigt-Reuss-Hill  method  (see 
discussion,  main  text ). 


A.  Seawater* 


Salinity 

(ppt) 

Density  Vp  Bulk  Modulus,  kh. 

(g/cc)  (in/ sec)  (dynes  cm2  X  1 0>  0) 

33.50 

1.0228  1 

1528.3  2.388955 

34.00 

1.0232  1 

1528  9  2.391766 

34.50 

1.0236  l 

1529.4  2.394266 

35.00 

1  024  ' 

1  $30.0  2.397082 

•Density  and  velocity  from  NAVOCLANO  SP-68  (1966);  computed  from:  p  V ^  =  k 

B.  Mineral  Species** 

Bulk  Modulus,  k 

Mineral 

(dynes  em-  X  10lf) ) 

Reteience 

C'aleite 

72.940 

Feselnick  (1962) 

Microcline 

5  1 .8  1  3 

Anderson  and  Safe  1 1965 ) 

Orthoclasc 

47.393 

Anderson  and  Nafc  (1965) 

Albite 

52.910 

Brace  <  1965). 

Anderson  and  Nafc  ( 1965 ) 

Labradorite 

66.667 

Birch  ( 1966) 

Quartz 

37.726 

Soga  (1968). 

McSkimin  cl  al.  ( 1965 ) 

Obsidian 

37.800 

Maiighnani  cl  al.  ( l%8 ) 

Hornblende 

84.175 

Brace  (1965) 

Biol  ilc 

4i .  i  52 

Brace  I 

Apatite 

91.743 

Birch  1 1966) 

Magnetite 

181.818 

Birch  ( 1966) 

Olivine 

126.582 

Birch  ( 1966) 

hnslatiie 

99.010 

Bnch  ( 1966) 

:  lyperstbene 

101.010 

. .  i  C  66) 

/vugite 

98  039 

Birch  ( 196' 

•’Computed  from  h  • 

-  l/,3  where  £  given,  references  are  u  •cent,  reuai!\  ^a.-ible  ..  fees  which 

frequ.'itlv  cite  other  sources;  for  recent, genera)  conipilutioni  sec  b‘r-:h  * !  ?n(> • .  -  Anderson 

and  Liebennunn  1 1968.1 


Computed  Aggieg.ilc 

Hulk  Moduli.  ks* 

Sediment 

Bulk  Modulus,  k 
(dynes/enr  X  10'® 

Ret  ci  once 

Medium  sand 

5  1 .206 

San  Diego  (nearshore) 

fine  gray  sand 

52.326 

San  Diego  (nearshore) 

Clayey  sill 

54.425 

San  Diego  Trough 

“Clay" 

50  -  -  - 

Skemplon  ( 1  %  1 ); 

used  for  deep-sea  sill- 
cla>s. 


•Some  typical  values. 
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APPENDIX  C: 

EQUATIONS  FOR  REGRESSION  LINES 
AND  CURVES  (ILLUSTRATED  DATA) 

Regression  linos  ami  cuivcs  wok  computed  t > >r  those  if  1  list r .t t col  sols  ol  ( v.i  )  ilala 
which  constitute'  the  best  indices  |  v)  to  oht.cn  desired  piopcrlies  (i  ).  Separate  equations 
ate  lisled,  wheie  appropiiale,  loi  each  of  the  three  general  environincnls,  as  follows:  eon- 
line  Dial  lei  nice  ( slid  I  and  slope  I.  (TV.  ahy  ssal  lull  (pelagic  I .  (II).  ahy  ssal  plain  I  tin  Indite), 

( P ) .  'Ihc  equations  aie  hosed  h\  Heine  ninnhois  to  the  lelated  scatter  diagrams  in  the 
in. tin  text .  The  Sland.nd  hnois  ol  I  stnnate.  n.  opposite  each  equation.  are  applicable 
mils  ncai  the  moan  of '  the  (,v,i  I  values,  and  acctu.tcs  ol  the  O')  value-,  given  (v ).  tails  oil 
away  fiont  lliis  icgion  (fiiiffiths,  1967,  p  448). 

It  is  iiupoitant  tlut  the  icgiessioii  ci|iiations  bo  used  only  between  the  limiting 
values  ol  the  index  piopcrty  (v  sallies),  as  noted  below  These  equations  aie  strictly 
empmeal  and  apply  only  to  the  (,v,i  )  data  points  involved. 

The  iiinitnig  values  ot  (  vl  in  tile  equations  below  aie  . 

1.  Porosity .  a,  peieent 
(T).  35  to  82  percent 

(.11)  and  (P).  30  to  40  percent 

2.  Density .  p,  g/cc 

(T),  1.25  to  2.10  g/cc 
(Hi.  I  25  to  1.50  g/cc 
(P).  1.15  to  1 .45  g/cc 

3.  Density  X  ( Velocity  )'.  p  I '  “  ,  dy  nes/ent2  X  1010 
(Hi,  2.4  to  3.4  dynes/cnr  X  I010 

(Pi.  2.7  to  3.4  dy  ncs/cm*  X  101  y 

4.  flay  siic  crams.  C\  percent 
(111,  45  to  80  percent 

frame  hulk  Modulus,  k  ,■  ( rlyncs/c m ~ )  vs.  I’qiqshy .  n  I  itaci  ion  )  l  I’lgine  2 ) 

(T I  For  sands:  k  dy  ncs/ein"  X  109 

logK(  =  2.71405  -  4.12135  ( .//) 

(T),  (H),  (P)  for  S lit -C i ;•>  s  s  dy  ucs/ui'r  X  I'd** 

logic.  =  3.73807  -  4.25571  (.in 

Bulk  Modulus,  K  ( dy  ncs/cm*  X  1 01  0  |  \s.  Pomsity  ,  n  (pctcenl )  (figures  7.  8 ) 

(  I )  c  =  22.73048  -  0.674982  (n)  +  0  0  0  8  265  (n)2  -  0.00003566  (n)3 
o  =  0.1379 

(11) tP)  x  =  4  05546-  0.1  14167  («) +  U.00047M/:)- 
o  =  0.014! 

Bulk  Modulus,  K  (dy  nes/cni2  X  1 0*  a  l  vs.  Density ,  p  U  cc )  ( figures  9.  1 0) 

(T)  K  =  -  22.16073  +  49  9361 16|p>  -  34. 144205  (p)2  +  8.1  19349(p)3 
u  =  0  1M2 


(II)  k  =  (>.09991  -  7.428005  (p)  +  3.445399  (p)2 

o  =  0.0302 

(l’|  k  =  3  42494  2.486980  (p)  +  1.596811  (p)2 

o  =  0.1525 

Hulk  Modulus, k  (dvuos/cnr  X  10,0)\s  LXmisiU  X  (Velocity  )2,pl'„2  (d> ncs/cm*  X 

- Tg^uo  l~ 

(II)  k  =  0.24037  +  0.891027  (pl^2) 

o  =  0.0528 

(P)  k  =  0.06033  +  0.74  23  09  (p>  p2) 

o  =  0.1446 

Rigidity  Modulus,  p  (dy  ncs, 'cur  X  IQl0)vs  Clay  Size,  C(pcn:ein)  (Figure  15) 

(H)  p  =  2.38685  -  0.105307(C)  +  O.OOlbl  1  (C)2  -  0  00000828  (O* 
u  =  0  0334 

Rigidity  Modulus,  p  (dviies/cm2  X  IQ1  u)  vs  Density  X  (Velocity  |- ,  p  l'„ -  (dynes/cm 

X  1 0 1  Q  )  (I  igure  1 0 1 

(P)  p  =  4.78616  -  5.154326  (pKp2j  *  1.797785  [pl'p1}* 

-  0.200143  (p»p2) 3 

o  =  0.0229 
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APPENDIX  D: 

ADDITIONAL  SEDIMENT  PROPERTIES 

1  ablcs  I  and  2  ttom  iVt  !(TP  143)  arc  included  in  this  appendix  as  a  conven¬ 
ience  lu  the  reader .  These  tables  furnish  additional  information  on  properties  of  the  sedi¬ 
ments  discussed  in  this  report  (TP  144)  Tables  D-l  and  D-2  can  be  corielated  with  tables 
I  and  2  in  the  main  text 
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Utbflr.'itOT)  values  2i  (  .  1  atmosphere 

Density:  s  iturated.  bulk  density .  porosily .  salt  tree;  ratio,  velocity  in  sediment/velocity  in  seawater  at  2V>(\  !  atmosphere,  and 
salinin  o'  sediment  pore-water. 

SI  standard  error  of  the  mean. 
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